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ABSTRACT 


The  stability  of  one-dimensional,  steady  detonations  to  periodic 

disturbances  transverse  to  the  flow  is  examined  in  the  limit  of  small 

wavelength,  i.e.,  2rt/e  -*  0,  It  is  found  that  any  unstable  disturbance 

grows  in  time  t  as  exp  (xt)  with  Im(T)/e  and  Re(T)  of  order  e°.  The 

asymptotic  criterion  for  stability  is  found  to  depend  largely  upon  the 

2  2 

steady-state  profile  of  c^  -  u  (where  c^  is  the  frozen  sound  speed  and 

u  is  the  mass  velocity  relative  to  the  von  Neumann  shock)  as  a  function 

2  2 

of  distance  behind  the  shock.  Detonations  for  which  c^  -  u  decreases 
monotonically  are  found  to  be  stable  (in  the  e  -*  »  limit),  but  stability 
in  cases  in  which  this  quantity  increases  either  monotonically  or  up  to 
a  maximum  (and  decreases  beyond)  is  determined  through  simple  integral 
functions  of  the  steady-flow  variables.  More  complicated  profiles  are 
not  treated  explicitly.  In  contrast  to  the  labor  involved  with  applica¬ 
tion  of  the  general  theory  of  detonation  stability,  the  current  asymptotic 
result  can  be  readily  applied  to  any  detonation,  irrespective  of  the 
number  of  chemical  reactions  which  occur,  provided  knowledge  of  the 
equation  of  state  and  reaction  kinetics  is  at  hand. 
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The  criterion  is  applied  to  an  idealized,  one-reaction  (A  -*  B) 

detonation.  Unstable  regimes  are  found  only  if  the  ratio  of  specific 

heats  y  is  less  than  2.  for  which  case  all  detonations  are  found  to  be 
'o 

unstable  for  sufficiently  large  values  of  the  activation  energy.  For  lew- 
heats  of  reaction,  it  is  found  that  instability,  although  not  reported  in 
an  earlier  n-umerical  application  of  the  general  theory  to  long  wavelengths, 
persists  for  all  (positive)  heats,  unless  the  activation  energy  is  also 
small.  On  the  other  hand,  for  more  realistic  heats  of  reaction,  the 
present  (asymptotic)  criterion  predicts  stability  for  small  activation 
energies,  where  long- wavelength, unstable  modes  were  previously  foimd. 
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GLOSSARY  OF  SYMBOLS 


This  glossary  contains  the  principal  symbols  used  in  the  main  text 
of  this  report,  hut  does  not  list  those  occurring  only  in  the  appendices. 
Symbols  which  are  not  included  are  components  of  vectors  and  matrices; 
these  symbols  are  subscripted  versions  of  the  principal  symbols,  with 
the  overhanging  arrow  (^)  for  a  vector  or  underscore  (*•)  for  a  matrix 
omitted.  In  the  following  list,  the  superscripts  and  subscripts  i  and  j 
are  understood  to  take  on  any  positive  integer  (or  zero)  value,  while 
V  denotes  any  signed  integer  or  zero.  All  other  subscripts  and  super¬ 
scripts  have  a  specialized  meaning. 


Symbol 

a^(z) 

A(z,  oo) 


A,(.)l 

V'- 

€)  i 

b^(z) 


Section 

Description  Introduced 

Polynomial  coefficients  in  the  solution  of  the  6 

second-order  related  equation. 

Polynomial  in  l/^>  with  coefficients  a^(z),  6 

entering  the  solution  of  the  second-order 
related  equation. 

Matrix  functions  of  the  steady  flow  which  enters  2 
the  ^  differential  equation. 

Integrals  over  the  reaction  zone  involving  2 

steady-flow  variables  and  6. 

Polynomial  coefficients  in  the  solution  of  the  6 

second-order  related  equation. 

Polynomial  in  i/cd  with  coefficients  b^(z), 

entering  the  solution  of  the  second-order 
related  equation. 


B(z,  to) 
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GLOSSARY  OF  SYMBOLS 
(Continued) 


c(x,  e) 


C 


C(x,  e) 


d(x) 

dg! 

e(x,  e) 


e, 

1 


e^(x,  1?) 


Matrix  function  of  the  steady  flow  which  enters  2 
the  ^  differential  equation. 

Frozen  sound  speed  in  the  steady  flow,  3 

Coefficients  of  admissible  asymptotic  solutions  4 
in  the  continuation  problem  for  0, 

Polynomial  coefficient  in  the  solution  of  the  6 

(n  +  2)-order  related  equation. 

Polynomial  in  l/e  with  coefficients  c^(x)  6 

entering  as  factor  in  p(x,  e) ,  the  solution 
of  the  (n  +  2) -order  related  equation. 

Contour  of  integration  for  the  definition  of  4 

a(e,  p). 

Matrix  of  n  +  2  linearly  independent  columns  6 

c^(x,  e). 

Derivative  of  c^t^  with  respect  to  x  in  the  steady  4 
o 

flow. 


Coefficients  of  admissible  asymptotic  solutions  4 
in  the  continuation  problem  for  0. 

Diagonal  matrix,  similar  to  0*  3 

Correction  term  in  the  second-order  related  6 

equation;  the  element  of  matrix  e). 

Unit  vector,  having  i^^  element  only,  3 

Leading  term  in  the  series  in  l/e.  3 


10 


E(x,  v>) 
f 


gi(x) 


h^(x,  0 


h^jCx.  0 


H(x) 

i*A 

i 


I(x) 


12(2) 


GLOSSARY  OF  SYMBOLS 
(Continued) 

Activation  energy  for  the  idealized  systems. 

A  transformed  version  of  ^^(x,  o) . 

Degree  of  overdrive  of  the  steady  idealized 
detonation. 

First  correction  term  in  the  series  in  l/e. 

Exponential  coefficient  of  l/e^  in  the  trans- 
formation  sequence  from  0  to  t. 

Functions  of  the  steady  flow  which  enter  the 
integrands  of  'b^(T,  e)  and  e) . 

Combination  of  the  heat  of  reaction  Q  and  the 
heat  capacity  ratio  7^  for  the  idealized 
system. 

Coefficient  of  e  in  the  exponential  of  the 
asymptotic  series  for 

b^(x,  C)  -  b^(x,  i). 

Vector  functions  of  the  steady  detonation, 
entering  into  V(t,  e). 

Regions  in  the  complex  z-plane. 

Transformed  version  of  £^(x,  u) . 

Generic  designation  for  the  integrand  for 
b^(T,  e)  and  b^CT,  e) . 

Integrand  of  the  expression  for 
idealized  system. 


.7 

3 

7 

3 

6 

2 

7 

3 

3 
2 

4 
6 

5 

7 
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GLOSSARY  OF  SYMBOLS 


(Continued) 


I.  Unit  matrix  of  order  i,  6 

^n+4*  ^ 

j(z)  Function  of  longitudinal  distance  entering  the  6 

solution  of  the  second-order  related  equation. 


e) 

Correction  term  in  the  (n  +  4)-order  related 
equation . 

6 

£ 

Pre-exponential  factor  in  the  Arrhenius  rate 

constant  for  the  idealized  system. 

7 

k. 

1 

Coefficient  of  in  k(u)) . 

6 

k(») 

Polynomial  in  l/cu  with  coefficients  k^  entering 
the  solution  of  the  second-order  related 

equation. 

6 

0 

0 

The  n  independent  part  of  kQ. 

6 

^01 

The  coefficient  of  0  in  kQ, 

6 

Coefficient  of  in  the  exponential  of  the 

asymptotic  series  for 

3 

kf(x,  t,  -o) 

Part  of  k^(x,  n)  which  remains  finite  at  the 

5 

turning  point. 


K 

^0 

K(z) 


Path  of  integration  from  x  =  0  to  ®  ,  5 


Quantities  independent  of  e  and  ,  6 

Factor  of  12(2)  for  idealized  system.  7 
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GLOSSARY  OF  SYMBOLS 
(Continued) 


^(x) 

\(x) 

VO  I 

LgCO  ) 

L(e,  t,x>) 


L 

HM& 


e) 

L(x,  e) 

m 


M 

f!(x,  €) 


n 


n 


(x) 


N.(x) 


Function  of  the  steady  flow;  its  value  behind  the 
shock  enters  L^  and 

Element  21  of  the  matrix  -  L. (x) ,  entering  the 
second-order  related  equation. 

Contributions  to  L(€,  Ob)  in  the  asymptotic 
limit. 

V(t,  e)/^  for  T  =  eC  +  n. 

The  Onsager  matrix  of  order  n. 

Coefficient  of  in  the  power  series  for  the 
coefficient  matrix  for  the  (J  ^  i) 

differential  equation. 

i-1  .. 

Polynomial  ^  e  ‘^L .  (x)  . 
d=0 

£W(X,  €). 

An  arbitrary  integer,  greater  than  0,  but  fixed 
throughout  the  discussion  in  Section  6. 

Generic  symbol  for  a  large  positive  number. 

Confluent  hypergeometric  functions. 

Lower  (n  +  2)-order  block  of  J(x,  e). 

The  integer  number  of  independent  chemical 
reactions. 

Function  of  the  steady  flow  entering  N^(x). 

Lower  (n  +  2)-order  block  of  ^^(x). 


5 

6 


2 

3 

6 


6 

6 

5 

6 
6 
1 

6 

6 
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GLOSSARY  OF  SYMBOLS 
(Continued) 


P(x) 

Po(>=) 

P]_(z) 

P^(z) 

P(x,  T,  e) 

q(x,  e) 


5(7o) 

Q 

Q(z,  oi) 

r(x) 

R  (x) 
R(x,  e) 

s(x,  0 
S(x) 


Pressure  in  the  steady  detonation. 

Leading  coefficients  in  the  l/e  expansion  of  the 
second-order  related  equation. 

Paths  in  the  complex  z-plane  used  to  establish 
the  admissibility  of  solutions. 

Coefficient  matrix  of  the  fundamental  9 
differential  equation. 

Power  series  in  l/e  which  enters  the  second-order 
related  equation. 

Coefficient  of  in  the  power  series  for  Q(z,  ui). 

Transitional  value  of  Q,,  determined  by  7^. 

Reduced  heat  of  reaction  for  the  idealized  system. 

Reduced  activation  energy  for  the  idealized  system 

Transformed  version  of  q(x,  e) . 

Transformation  matrix,  which  is  part  of  ^(x,  €) . 

Chemical  rate  in  the  steady  detonation. 

Coefficient  of  in  the  power  series  for  ^(x,  e) 

Correction  term  in  the  (n  +  4) -order  related 
equation  for  cp(x,  u,  e) . 

2  2 

Squa.re  root  function  of  c^q  and  ^  . 

Specific  entropy  in  the  steady  flow. 

Sectors  of  the  complex  z-plane. 


GLOSSARY  OF  SYMBOLS 


(Continued) 


Transformation  matrix,  entering  the  sequence 

from  0  to 

6 

S^^^\x,  e) 

The  product  S_  •  S,  • • •  S . . 

^  <~>0  —I  '*»a. 

6 

£.(x,  e) 

6 

t 

Time  . 

2 

T(x) 

Temperature  in  the  steady  flow. 

5 

T(x,  0 

Transformation  matrix  which  diagonalizes  ^(x,  0  • 

5 

u(x) 

X  component  of  velocity  in  the  steady  detonation, 

5 

u(z,  cd) 

L  (z, 
v''  ■ 

u(x) 


U^(z,  ($)) 


U^^\x,  e) 


in  the  frame  of  reference  of  the  steady  shock. 

Solution  of  the  second-order  related  equation. 

Mass  velocity  in  the  steady  flow  with  components 
u(x)  and  0, 

Coefficient  of  e“^  in  the  power  series  for 
£^^\x,  e). 

Portion  of  coefficient  matrix  of  the 

differential  equation  not  in  the  block- 
diagonal  form  of  L^^\x,  e)  , 


6 

2 


U.(x) 

6 

e) 

--(m)  /  N 

U)  ^(x,  e). 

6 

v(x) 

Specific  volume  in  the  steady  flow. 

6 

&) 

Contribution  to  oS)  in  the  second-order 

related  equation. 

6 
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GLOSSARY  OF  SYMBOLS 
(Continued) 


V(t,  e)  Fundamental  function  whose  roots  in  the  T-plane 

determine  stability. 

W(t,  e)  Contribution  to  the  shock  distortion  from  the 

initial  perturbations . 

W^[^,  k(cD)]  A  solution  of  the  confluent  hypergeometric 

differential  equation. 


1 

2 

6 


X* 


X 

o 

X^(x) 

y 

y(x,  €) 
e) 

y(x,  e) 
y^Cxj  e) 


Whittaker  functions .  6 

Variant  of  W  x(^) •  6 

kj  4 

Longitudinal  coordinate  in  the  frame  of  reference  2 
of  the  steady  shock. 

Value  of  X  at  a  turning  point,  5 

The  greater  of  a  pair  of  turning  points.  4 

Pa,rticula.r  value  of  x.  3 

Function  of  the  steady  flow  entering  0((e,  'o).  6 

Spatial  coordinate  transverse  to  the  steady  flow,  2 

Component  1  of  y(x,  e)  or  y^(x,  e) ;  solutions  of  6 

the  second-order  related  equation. 

Solutions  of  the  (n  +  4) -order  related  equation.  6 

Spatial  variable  derived  from  x.  Except  in  4 


Section  6,  it  denotes  the  analytic  continuation 
of  X  to  the  complex  plane.  In  Section  6,  it 
represents  x  -  x*. 


l6 


<ts] 


GLOSSARY  OF  SYMBOLS 
(Continued) 


Variable  of  integration,  replacing  x,  used  in  5 

several  contexts, 
z(0  Variable  of  Integration  for 

idealized  system, 

«(€,  Coefficient  of  the  Hg-'type  contribution  to  5 

for  turning  points  associated  with  profiles 
M  and  I, 

Ci:(e,  13)  Function  having  the  properties  of  0!(6,  "o) ,  but  6 

for  the  turning  point  at  the  maximum  of 
profile  M. 

p(e,  x>)  in  0!(€,  p),  '5 

Sl(^i)  Coefficient  of  -i e  in  p(e,  p) ,  5 

P^C^i)  Contribution  to  p(e,  h)  independent  of  e  and  p,  5 

Pj(^i)  Coefficient  of-PinP(€,^,  p),  5 

Ratio  of  specific  heats  for  the  idealized  system,  7 
r(a)  Gamma  fimction,  6 

6  Either  i  or  1,  depending  on  whether  the  steady  6 


detonation  has  a  maximum  or  minimum  at  the 
extremum  of 

A(x,  Cj 

V") 

A2(^) 


Factors  in  the  asymptotic  expressions  for  the 
integrands  of  b^(T,  e)  and  b2(T,  e) , 

Wave  number  of  the  disturbance  (times  2-n')  in  the 
transverse  spatial  coordinate  y. 


IT 


GLOSSARY  OF  SYMBOLS 


(Continued) 


^1 

Tl(x) 

e(x,  T,  e)  ) 
e(x,  o,  e)  j 


Leading  term  in  the  assumed  dependence  of 
T  on  e;  t  =  +  o. 

Im(0. 

Sonic  parameter  in  the  steady  flow. 

Bounded  solution,  as  x  of  the 

fundamental  differential  equation. 

Solution  of  the  fundamental  differential 

equation,  irrespective  of  the  boundedness 
condition. 


|(v) 


Kix) 


g)  Series  expression  for  analytically 

continued  axound  a  turning  point  to 
sector  S^. 

Mach  number  in  the  steady  flow. 


K  Shock  Mach  number  for  idealized  detonations 

marking  the  transition  from  profile 
M  to  profile  I  behavior. 


2 

2 

3 

2 

3 

4 

3 

7 


\(x) 


k(x) 

u) 


Progress  variable  in  the  steady  flow  for  the 
idealized  system. 


Progress  variable  in  the  steady  flow. 


Coefficient  matrix  for  the 
equation. 


differential 


Coefficient  of  in  the  power  series 

expansion  of  the  eigenvalues  of 
-  P'(“5  e)  in  l/e. 


7 

2 

6 

4 
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GLOSSARY  OF  SYMBOLS 
(Continued) 


5) 

V 

V 

») 

?(>.) 

e) 

n(x,  C,  ^5  e) 

rtj_(x,  C,  ^5  e) 
n(z) 

P(x,  e) 

a(v) 

a(x) 

Z 

V 

Z(z,  &) 


Eigenvalue  of  ^(x,  ^)  . 

^^j_(x,  0  -  0* 

Compression  ratio  v^/’v  fo^  "th®  steady  shock. 

Coefficient  of  e"^  in  the  power  series 
expansion  of  the  eigenvector  of 
-  P'(“,  +  ^5  e). 

Monotone  function  of  the  progress  variable 
for  the  idealized  system. 

Value  of  ^\)  at  chemical  equilibrium;  ^O). 

Value  of  C(^)  at  the  turning  point. 

Fourier-Laplace  transform  of  the  shock 
distortion. 

T(x,  0"^  “  9(x,  e). 

T(x,  0"^  •  9.(x,  o,  e). 

Factor  in  the  solutions  of  the  second-order 
related  equation,  u^(z,  e) . 

Solution  of  the  (n  +  2)-order  related  equation. 

(-1^. 

Thermodynamic  function  in  the  steady  flow. 

Sector  in  the  complex  f’-plane. 

Combination  of  coefficients  in  the  second-order 
related  equation. 


5 

4 

5 

4 

7 

7 

7 

2 

5 

5 

6 

6 

6 

3 

6 

6 
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GLOSSARY  OF  SYMBOLS 
(Continued) 


T 

•D 

0. 

1 

T(z,  cS) 

9 

9 

<t(z) 

t,  e) 

\|/(t,  €) 

^(y?  t) 
t(z,  &) 

e) 

t(x,  e) 

?^(x,  e) 


Laplace  time  transform  parameter. 

Part  of  T  assumed  independent  of  e, 

Im('o) . 

Factor  in  the  solutions  of  the  second-order 
related  equation  u^(z,  &)  . 

arg( z  -  X*) . 

arg(x*  -  z) . 

Po(^)^. 

Integrals  of  l-i^(x,  and  ^2(^5  ’  entering 

the  asymptotic  expressions  for  b^(T,  e) 
and  bg  ( T ,  e)  , 

Solution  of  the  (n+ 4) -order  related  equation. 
Leading  term  in  e  of  -  ^) » 

Part  of-P'(x,  ^e+u,  e)  independent  of  e. 
Fourier  transform  of  t)  • 

Distortion  of  the  shock* 

Stretched  coordinate  in  the  second-order 
related  equation* 

Transform  of  0(x^  o,  e)  by  e) , 

Independent  solutions  of  the  differential 
equation . 


1 

2 

6 

6 

4 

4 
6 

5 

6 
2 
2 
2 
2 
6 

6 

6 

6 
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GLOSSARY  OF  SYMBOLS 


(Continued) 


Parameter  for  the  idealized  system  determined 
by  detonation  velocity  and  7^. 

Large  parameter,  6e, 

Parameter  for  the  idealized  system  determined 
.  by  detonation  velocity  and  7^0 

Interval  of  the  real  x  axis  on  which  T(x, 
is  differentiable. 


Subscripts 

Evaluation  in  the  steady  flow,  behind  the  shock. 

Evaluation  in  the  steady  flow,  in  front  of  the 
shock. 

Evaluation  at  the  point  in  the  steady  flow  where 

2  .  . 

c  n  IS  maximum, 
o 

Evaluation  for  the  Chapman-Jouguet  detonation 
velocity. 


1.  imRODUCTION 

The  phenomenon  of  detonation  has  been  the  subject  of  continuing 
theoretical  analysis  ever  since  the  pioneering  work  in  the  nineteenth 
century.  As  a  result  'Of  a  certain  amount  of  agreement  between  theory 
and  experiment,  this  effort  resulted  in  the  general  acceptance  of  the 
so-called  von  Neumann-Zeldovich-Doering  theory.  The  fundamental  notion 
in  this  theory  is  that  the  phenomenon  consists  of  a  shock  propagating 
into  the  quiescent  material,  with  a  resultant  initiation  of  exothermic 
chemical  reaction  behind  this  front,  A  quantitative  description  of  the 
structure  of  the  detonation  wave  is  presumably  complicated  in  its  early 
stages,  but  the  theory  directs  its  attention  to  the  situation  after  all 
transients  arising  from  the  initiation  of  the  detonation  have  vanished. 

It  is  supposed,  then,  that  the  detonation  ultimately  becomes  a  steady, 
one-dimensional  flow  when  viewed  in  the  frame  of  reference  of  the  shock. 
The  investigation  of  the  consequences  of  these  assumptions,  with 
particular  attention  to  the  Chapman-Jouguet  hypothesis  for  the  detemnina- 
tion  of  the  velocity  of  the  shock  front,  has  been  a  principal  concern  of 
theoretical  endeavors  in  this  field. 

Beginning  approximately  in  I959  the  experimental  evidence  was  in¬ 
creasingly  seen  to  be  not  entirely  compatible  with  the  above  picture.  In 
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particular,  numerous  instances  of  three-dimensional  structure  of  an 

1  2 

apparently  non-transient  nature  ha.ve  been  reported.  ’ 

From  a  theoretical  point  of  view,  the  possibility,  that  the 
supposed  one -dimensional,  steady  detonations  do  not  actually  occur,  is 

3 

subject  to  analytic  investigation.  The  author  has  attempted  to  deter¬ 
mine  the  so-called  hydrodynamic  stability  of  steady  detonations;  that  is, 
whether  or  not  small  disturbances  imposed  upon  the  steady  flow  at  some 
instant  will  die  out  with  time.  Should  any  infinitesimal  disturbance 
grow,  the  detonation  is  said  to  be  hydrodynamically  unstable. 

The  key  to  the  hydrodynamic -stability  theory  is  that  the  equations 
of  reactive  hydrodynamics  are  linearized  in  the  deviations  of  the  per¬ 
turbed  flow  from  the  steady  flow.  By  virtue  of  this,  the  dependence  of 
the  perturbations  on  the  coordinates  transverse  to  the  steady  flow  is 
removed  by  Fourier  decomposition  and  one  is  left  with  the  stability 
problem  for  a  given  transverse  wave  number  e/2jt. 

Now  the  general  theory  of  detonation  stability^  can  be  carried 
to  the  definition  of  a  criterion  based  upon  the  roots  in  the  complex 
T-plane  of  a  certain  function  V  of  t,  g,  and  the  steady  flow.  The 
determination  of  V(t,  g)  for  a  given  steady  flow  is  by  no  means  simple, 
requiring  the  solution  of  a  system  of  ordinary  differential  equations  of 
order  n  +  4  (for  detonations  involving  n  independent  chemical  reactions) 
whose  coefficients  depend  upon  the  steady-detonation  structure  in  the 
reaction  zone  as  well  as  t  and  g.  Ordinarily,  then,  the  stability 
problem  is  expected  to  require  numerical  solution  with  results,  therefore. 
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limited  to  a  finite  range  of  wa.ve  numbers.  A  recent  calculation  for  an 
ideal-gas,  unimolecular-reaction  system  was  able  to  establish  instability 
of  certain  detonations  over  finite  ranges  of  e.  However,  any  assertion 
of  the  stability  of  a  given  detonation  was  necessarily  qualified  to 
include  only  those  values  of  e  which  were  numerically  accessible  to  the 
calculation.  Thus,  we  are  led  to  investigate  the  stability  criterion  in 
the  limit  of  small  wavelength,  i.e.,  large  e. 

Since  the  present  theory  neglects  the  transport  effects  of 
diffusion,  viscosity,  and  heat  conduction,  the  unit  of  length  in  the 
problem  is  determined  by  the  gradients  in  the  steady  flow.  Hence,  the 

analysis  of  the  e  “  limit  is  expected  to  apply  to  wavelengths  which 

t 

are  short  relative  to  the  re action- zone  length.  On  the  other  hand, 
transport  effects  are  expected  to  become  important,  at  least  for 
finite  disturbance  magnitudes,  for  wavelengths  comparable  to  a  mean 
free  path.  Although  the  validity  of  our  theory  cannot  be  ascertained 
without  recourse  to  the  transport  equations,  it  is  not  unreasonable  to 
suppose  the  restriction  of  our  present  considerations  to  wavelengths 
which  are  long  relative  to  a  mean  free  path  and  short  relative  to  the 
reaction-zone  length  to  be  physically  meaningful. 

The  mathematical  considerations  which  are  attendant  upon  the 
investigation  of  the  e  “  limit  are  those  of  the  so-called  parameter 

The  reaction  zone  is  formally  of  infinite  extent.  The  phrase  "reaction- 
zone  length"  refers,  of  course,  to  any  convenient  measure  of  the  region 
over  which  state  variables  change  substantially,  e.g,,  the  distance  to 
half -completion  of  some  reaction. 
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proTjlem  for  linear  differential  equations.  This  branch  of  analysis  is 


essentially  complete  for  systems  of  differential  equations  of  any  order, 
provided  there  are  no  so-called  turning  points,  but  is  limited  to  certain 
second-order  equations  and  specific  high-order  equations  when  a  turning 
point  is  involved.  For  large  systems  of  equations  there  does  not  seem 
to  be  any  appreciable  analysis  in  the  literature  for  the  turning-point 
problem. 

Turning  points  of  rather  simple  types  play  a  central  role  in  the 
analysis  of  the  small-wavelength  limit,  so  that  we  require  that  the  small, 
but  troublesome,  gap  associated  with  the  (n  +  4)-order  system  of  equations 
first  be  bridged.  A  large  portion  of  this  report  is  involved  with  pre¬ 
cisely  this  problem. 

After  establishing  asymptotic  expressions  for  the  desired  function 

/  ^  4 

V(t,  e),  the  stability  of  the  idealized  system  previously  investigated 
is  examined  in  the  small-wavelength  limit. 

2,  resume'  of  the  geheral  theory 

The  time  evolution  of  an  initial  perturbation  from  a  steady 
detonation,  moving  along  the  x  axis,  is  reflected  by  the  secular  beha,vior 
of  the  distortion  if(y,  t)  of  the  leading  shock.  If  x  is  the  coordinate 
in  the  direction  of  the  steady-state  flow  and  the  origin  (x  =  O)  is 
located  in  the  (distorted)  shock,  while  y  is  the  coordinate  transverse 
to  the  steady  wave,  then  -  ?(y,  t)  is  the  x  coordinate  of  the 
1 - ; - 

For  brevity,  we  use  only  one  transverse  coordinate.  The  stability 
criterion  is  identical  to  that  for  three  dimensions. 5 
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unperturbed  shock  at  time  t.  The  general  theory  of  detonation  stability 
has  shown  that  for  infinitesimal  disturbances  the  Fourier  transform 
of  ^  on  the  transverse  coordinate  has  its  Laplace  time  transform 


>|-00 

^(■tj  =  f  ip(y,  t)  dy  (2.l) 

«  00 

given  by 

|(t,  e)  =f  e"^iKt,  e)  dt 

‘'o 

=  W(t,  e)/v(T,  e) »  (2.2) 

The  function  W(t,  e)  is  determined  by  the  e-Fourier  component  of  the 
initial  perturbations  in  the  state  variables  [say,  specific  volume  v, 
mass  velocity  u  (a  2-vector),  specific  entropy  S,  and  chemical  composi¬ 
tion  \  (an  n-vector)]  throughout  all  of  space,  and  has  been  shown  to’ 
contain  no  singularities  in  the  right  half  x-plane. 

The  function  V  depends  upon  the  steady  flow  but  not  on  the 
initial  perturbations.  It  is  given  explicity  by 

V(t,  e)  =  rb^(T,  e)  +  ieb2(T,  e)  -  0(O,  t,  e)  *  (Th^  + 

/“  -*  1  -* 

e(x,  T,  e)  •  ^  (x)  *  g^(x)dx  (2.3) 

0 

00 

The  vector  function  0(x,  x,  e)  is  the  solution  of  the  differential  equation 
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(2.4) 


^  -  P'(x,  T,  e)  •  e 

P(x,  T,  e)  =  -  +  i€Ay(x)  +  B(x)] 


which  is  bounded  in  the  reaction  zone  for  Re(T)  >  0;  i.e., 


0(x,  T,  e)  =  0(l)  as  X  “  for  Re(T)  >  0.  (2.5) 

The  prime  denotes  the  transpose  matrix.  The  vectors  0,  g^,  gy>  h^,  and 
h  are  (n  +  4) -vectors,  one  component  for  each  state  va.riable,  while 

y 

A  (x),  A  (x),  B(x),  and  I  are  (n  +  4)-order  square  matrices;  ^ is  the 
unit  matrix  and  the  remainder  of  these  (Quantities  (given  in  Ref*  3)  are 
completely  determined  by  the  steady  detonation,  with  h^  and  h^  being 
determined  solely  by  the  von  Neumann  shock  itself,  i.e,,  by  the  equation 
of  state  of  the  unreacted  explosive  and  the  steady  detonation  velocity. 
The  occurrence  of  instability,  in  the  sense  of  exponentially 
growing  contributions  to  'k(t,  e) ,  is  seen  to  be  occasioned  by  roots  of 
V(t,  e)  lying  in  the  right  half  x-plane.  The  problem  of  detonation 
stability  resolves  into  determining  whether  or  not  V  has  such  roots. 

The  only  general  result  in  this  regard  is  the  fact  that  for  [t  |  large 
(at  fixed  values  of  the  detonation  velocity  and  e)  with  Re(T)  S;  0,  V  is 

p 

of  order  |t  |  and  has  no  roots  at  infinity. 

Now  in  obtaining  V(t,  e)  asymptotically  in  e,  it  is  of  importance 
that  the  magnitude  of  t  be  permitted  to  depend  on  g,  the  dependence 
being  such  as  to  reproduce  the  asymptotic  dependence  of  the  roots  of  V  on 
G,  For  this  purpose,  we  let 
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(2.6) 


T  =  +  V) 

whence  Eq,  (2.5)  3,nd  (2,4)  becoine 

V(t,  e)  =  eL(e,  C, 

L(e,  +  1^2  -  e(0>  U  e)  *  ^ 

+  €“^u  [b^  -  e(0,  u,  e)  •  h^] 

g  =  [€^(X,  0  +  O)]  •  5 

Jo=  •  (5i*  ^iy”' 

ii=  i:£^  •  ("J.-^s.)]' 

?  =  0(l)  as  X  for  Re(T)  >  0  .  (2.Tc) 

The  significance  of  the  substitution,  Eq,  (2.6),  is  that  the 
leading  term  of  9  (as  e  -*  ®)  depends,  as  we  shall  see,  upon  both  and 
0  ,  and  any  unstable  foots  of  the  leading  term  of  V(t,  c)  occur  on  the 
axis  Re(^)  =  0,  Thus,  the  sign  of  the  real  part  of  O  will  determine 
whether  such  roots  are  stable, 

5.  FORMAL  SOLUTIONS 

In  this  section,  we  obtain,  following  the  theory  of  the  parameter 
problem,^  the  leading  terms  in  the  e-asymptotic  solutions  of  Eq,  (2,7b), 
The  validity  of  these  expressions  in  the  asymptotic  representation  of 
^  will  be  the  concern  of  Section  4, 
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Series  solutions  of  Eq.  (2.7b)  are  obtained  from  diagonalization 
of  the  leading  matrix  If  there  exists  a  similarity  transform 

T(x,  differentiable  on  an  x  interval  such  that 

(the  second  equality  denoting  that  ^  is  diagonal  and  defining  its 
elements),  then  the  transformation  of  9  by  ^  ^  yields  the  system  of 
equations, 

—I 

rt  =  T  •  9 


£1 

dx 


(eD  +  E) 


jr 


(5.2) 


^(x,  n)  =  T(x,  0"^  •  *  T(x,  0 


dT(x,  O/dx]. 


5  •  ' 

Follo'wing  the  scheme  given  by  Friedrichs ,  we  substitute 

^i  ^  ^)  +  u)  +  ...]  exp  [eh^(x,  i)  +  k^(x,  u)] 

(5.5) 

in  which  subscript  i  distinguishes  among  the  solutions  which  will  arise. 

The  unknowns  h. ,  e. ,  k. ,  f .  ...  are  found  by  equating  powers  of  e  in 
X  1  1  1 

Eq.  (5.2);  one  obtains 


4*  ^ 

The  factor  exp  (k.  )  is  not  separated  from  e^  and  f^  by  Friedrichs  j  but 
this  is  not  a  consequential  difference. 
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(5.4) 


X 

but  the  determination  of  the  remaining  quantities  depends  largely  upon 

—4 

■whether  or  not  the  eigenvalues,  are  distinct.  The  solution  is 
said  to  be  of  the  type. 

The  eigenvalue-eigenvector  problem  for  ^  can  be  solved  explicitly 
and,  in  fact,  is  similar  to  that  which  occurs  in  the  problem  of  shock 
stability.^  In  Appendix  A,  we  compute  ^  and  find  its  eigenvalues  to  be 

p^  =  -  /c(/c^  +  s)/tiu 

p  =  -  K(Kt,  -  s)/tiu 

(3.5) 

Pj  =  Pi^  =  . . .  =  =  ^/n 

s=[^^  +  c2tiF. 

The  notation  of  Eq.  (3.5)  is  as  follows:  t](x)  is  the  sonic  parameter 

n  2/2 

Ti  =  1  -  u  /c^ 

in  which  c^(x)  is  the  frozen  sound  speed  and  u(x)  is  the  mass  velocity  in 
the  steady  frame  of  reference;  k(x)  is  the  local  Mach  n-umber  ^a/c^.  The 
square  root  s(x,  ^)  is  to  be  ■understood  as  a  function  of  ^  ■vd.th  parametric 
dependence  of  x,  whence  it  suffices  to  specify  that  s  is  the  positive 
branch  of  the  function  ■with  a  branch  cut  along  the  imaginary  ^  axis  between 

±  ic  The  ina,trix  of  right  eigenvectors  T  and  its  inverse  are  also 

o' 

given  in  Appendix  A,  along  ■with  certain  elements  of  ^  which  will  be 
required. 
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Now  in  the  subsonic  reaction  zone,  we  see  from  Eq,  (5*5)  that  only 
IJ.^  has,  for  Re(^)  >0,  negative  real  part.  It  is  not  surprising  that 
will  be  of  primary  importance  in  finding  0  since  the  latter  must 
satisfy  the  boundary  condition  Eq,  (2.7c),  The  determination  of  the 
unknowns  in  is  straightforward; 

^^(x,  P)  = 

dk^(x,  ^)/<ax  =  E^^(x,  p)  (5.6) 


k^(x,  i,  ^)=f 


where  e  .  denotes  the  unit  vector  having  only  its  component  non-zero 


0 

The  lower  limit  of  integration  x^  need  only  be  selected  on  the  interval 
in  order  to  ensure  existence  of  the  integral.  The  next-order  term  is 


found  to  be 


P) 


0=2  o 


E2i/(M2  - 


(5.7) 


where  x  is  not  necessarily  the  same  as  in  Eq,  (5.6)  but  does  lie  on 
o 

1J<.  Higher  terms  are  readily  written  down  for  this  series. 


r 
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Now  it  will  be  convenient  to  denote  the  solutions  of  Eq,  (2.7b), 

irrespective  of  the  boradary  condition  Eq.  (2.7c),  by  subscripted 

symbols  9^  as  determined  by  the  inverse  transform  from  As  we  noted, 

then  9^  is  bounded  in  x  for  Re(0  >  For  Re(0  =  0>  1^^  has  negative 

1. 

real  part  for  x  large  only  if  |lm(0  |  ^'^o^^^x-co*  pure  imaginary  ^ 

of  greater  magnitude,  it  is  seen  that  9^  is  bounded  if  E^^  has  a  non¬ 
vanishing,  negative,  real  part.  From  Appendix  A,  Eq,  (A.16),  we  find, 

-4  .7 

using  linearization  of  the  reaction  rate  r  near  equilibrium. 


V”’  =  I  -  [ 


20  +  c  a  •  L 

O 


=]| 


(5.8) 


where  L  is  the  symmetric,  positive-definite,  Onsager  matrix  (order  n) 


and  a  is  a  thermodynamic  function,  related  to  the  heat  of  reaction. 


7 


It  follows  that,  for  Re(0  =  0,  U  |  >  (c  'ri^)„  ,  and  Re(o)  ^  0,  we  have 

O  X — 00 


Re[E^^(“>,  -o)]  <  0.  (5.9) 

Therefore  5  0^  is  hoimded  as  x  “  for  all  Re(T)  ^  0  [  except  ^  of  course , 

1 

for  t  =  ±  i(c  .  for  -which  the  interval  of  definition  must  he 

hounded  a.va.y  from  x  =  oo]* 

Although  0^  has  the  desired  boundedness  property^  it  will  he  seen 
that  02  is  also  required  explicitly*  The  evaluation  of  the  latter  is 
evident  from  our  expressions  for  9^.  Thus 
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(5.10) 


o 


The  other  solutions  0^,  6^}  •«•>  ®n+4’  ^5  need  not 

concern  us  other  than  the  fact  that  their  series  exist;  the  prescription 
given  by  Friedrichs^  applies  to  this  degenerate  case  as  well,  although 
explicit  determination  of  rt^,  ...  in  terms  of  elementa,ry  functions  is 
not  generally  possible. 

The  region  of  definition  of  the  9^  depends  upon  the  existence  and 

differentiability  of  the  transfom  ^(x,  0.  As  is  seen  from  Appendix  A, 

the  intervals  can  be  terminated  only  at  a  point  x*  at  which  s(x*,  C) 

vanishes  or  where  u(x*)  =  It  is  convenient  then  to  classify  the 

values  of  C  according  to  the  properties  of  the  matrix 

1  '  1 

Class  III:  Re(0  =  0  and  min  (c^n^)  ^  U  I  ^ 

X  .  X 

Class  II:  lm(0  =  0  and  min  (u)  ^  ^  max  (u) ,  _  (3.1l) 

X  X 

Class  I:  All  other  C?  with  Re(0  ^  0. 

For  ^  in  Class  I,  a  single  interval  extends  from  x  =  0  to  In  Class 
III  and  Class  II,  the  x  axis  is  broken  up  into  two  or  more  open  intervals j 
and  we  denote  these  •••»  beginning  at  the  right.  It  is  of  interest 

to  note  that  in  Class  II,  the  solutions  Gg  and  0^  become  "identical"  at 
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X*  in  the  sense  that  the  eigenvectors  of  ^  (see  Appendix  A)  correspond¬ 
ing  to  become  identical,  as,  of  course,  do  the  eigenvalues 

themselves.  Similarly  in  Class  III,  0^  and  become  identical  at  x*. 

4.  ASYMPTOTIC  EVAIUATION  OF  0 

The  desired  solution  0  is  defined  by  the  boundedness  condition, 

Eq.  (2,7c).  For  la.rge  x,  the  behavior  of  the  solutions  of  the  differ- 
ential  equation  (2.7b)  was  shown  previously  to  be  exponential  growth 
or  decay,  the  coefficients  of  x  in  the  exponentials  being  just  the 
eigenvalues  of  the  coefficient  matrix,  evaluated  at  x  =  “.  Moreover, 
for  Re(T)  positive,  it  was  shown  that  a  single  eigenvalue  has  negative 
real  part. 

The  formal  solutions  Eq,  (3.5) j  derived  in  the  last  section, 
behave  in  the  same  fashion  for  large  x,  for  the  h^  and  integrals 
become  linear  in  x  near  the  region  of  chemical  equilibrium.  Now  for 

Q 

large  €,  the  eigenvalue  problem  at  x  =  ”  has  a.  perturbation  solution 
through  a  power  series  in  e“^, 

[g  ^(»,  0  +  ^^(co,  v)  -  (ed^°^  +  +  ..•)ll 

•  +...]=  0 

which  exists  provided  that  t)  has  distinct  eigenvectors.  Then, 

with  the  exception  of  the  case  of  degenerate  eigenvectors  of 
[viz,,  f  =  ±  i(c  ri^)  and  t  =  (u)  which  we  do  not  treat  directly], 

’  '  O  X=oo  X=ra 
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it  Is  seen  tha,t  the  x  ->  «>  solutions  agree  with  the  expressions  of  the 
last  section  when  carried  to  large  x.  The  e  -*  oo  limit  and  the  x  “ 
limit  are,  therefore,  interchangeable 5  and  the  solution  6  can  be  investi¬ 
gated  directly  through  the  x  -  “  behavior  of  the  e-asymptotic  solutions. 

The  validity  of  the  expressions,  Eq.  (3.5))  as  asymptotic  solutions 
of  the  differential  equations  (2.7b)  depends,  according  to  the  theory 
of  the  parameter  problem,^  only  on  the  behavior  of  the  functions  h^(x,  . 

It  will  be  convenient  in  what  follows  to  regard  the  formal  solutions  as 
defined  not  only  on  the  axis  of  real  values  of  x  but  also  for  complex 
values  z.  For  e  real,  then,  the  region  of  validity  (the  admissible  region) 

ijTi 

K  of  the  z-plane  for  the  i  solution  depends  upon  the  existence  of  paths 
i 

P.(z)  in  K.  for  every  z  such  that  P.  is  of  bounded  contour  length  and 
J  1  d 

such  that  h.  -  h.  is,  in  real  part,  non-increasing  along  the  path,  i.e., 
d  1 


Re[h..(z')  -  h..(z)]  ^  0 


h..  =  h.  -  h. 


di  d  1 


(it.i) 


for  all  z'  on  P.,  from  an  origin  z'.  to  the  point  in  question  z.  The 
t]  3 

existence  of  such  paths  for  every  J  is  sufficient  to  prove  the 
a,(imissibility  of  solution  i  in  JC  • 

In  order  to  consider  condition  (4.l)  in  connection  with  G  ,  we 


write 


h,n(z)  =  li),n(z)  =  ...  =  h 


‘31V -  ...  -  %+4,iv“/  -31 


1^,1  (zOiz' 


(4.2) 


p.  .  =  n.  -  p . 

Id  1  d 
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We  consider  and  along  the  real  axis  for  the  present,  restricting 
attention  to  an  interval  for  which  9^,  62 f  •••  a'^e  defined.  As  long  as 
ReCUg^)  and  Re(p^^)  do  not  change  sign,  ReCh^^)  and  Re(h^^)  are  monotone 
in  X,  If  Re(h^^)  is  monotone  non-increasing,  then  a  suitable  path 
originates  at  the  left  and  extends  to  increasing  values  of  x.  If  Re(h^^) 


is  monotone  non-decreasing,  then  a  suitable  path  originates  at  a 
large  value  of  x.  and  extends  to  the  left.  In  either  case,  an  admissible 
region  3<^  can  be  found,  which  includes  that  portion  of  the  real  axis  on 
which  neither  Re(M-2-j_)  nor  Re(lJ.^^)  vanishes.  In  view  of  the  interchange- 
ability  of  the  e  and  x  limits ,  it  is  clear  that  such  a  region  can  be 
extended  to  infinite  x. 

From  Eq,  (5.5)?  we  compute 


^21  =  2/cs/tiu 

=  (t  +  Ks)/t\\x,  (4.5) 

Now  for  any  it  is  seen  that  Re(lJ.2^)  and  Re(Pj^)  do  not  change  sign  in 
any  interval  of  definition  of  0^,  62?  .•*>  fon  the  points  where  s 
vanishes  are  just  the  endpoints  of  the  intervals  for  Class  III  values  of 

Similarly,  the  validity  of  is  proved  by  consideration  of 
=  -  2/cs/tiu 

IJ.J2  =  (^  -  '^s)/tiu  (4.4) 

whose  real  parts  either  vanish  identically  or  do  not  vanish  at  all 
(depending  on  ^)  on  the  bb  intervals. 
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For  ^  in  Class  I,  it  follows  that  the  admissible  region 
includes  the  entire  x  axis  so  that 

^  ~  0^  for  all  X,  (4,5) 

For  Class  II  and  Class  III,  however,  we  have 

0  ~  0^  for  X  on  (4,6) 

but  are  left  with  the  problem  of  the  continuation  of  0  to  the  left  of  the 
so-called  turning  points  x*. 

There  are  essentially  two  techniques  used  in  solving  turning-point 
problems  and  it  seems  appropriate  to  pause  to  mention  some  general 
characteristics  of  these.  The  first  involves  an  attempt  to  proceed 
around  the  turning  point  by  analytic  continuation  into  the  complex  z 
(independent  variable)  plane.  Ordinarily  this  process  is  terminated  at 
certain  Stokes  lines  across  which  the  asymptotic  form  of  the  solutions 
to  the  differential  equations  changes  abruptly.  Usually  this  approach 
involves  solution  of  the  "connection  problem;"  that  is,  it  involves 
the  determination  of  the  solutions  at  the  turning  point  itself  by 
means  of  an  appropriate  "stretching  transformation,"'^  This  first 
approach,  although  extensively  investigated,  seemingly  has  not  been 
developed  to  any  degree  of  generality,  with  proofs  of  its  validity 
being  not  entirely  satisfactory. 

A  second  approach  is  that  of  Langer  and  McKelvey^~^  which 
concentrates  at  the  outset  on  the  nature  of  the  equations  at  the  turning 
point,  in  much  the  manner  of  the  solution  of  the  connection  problem,  but 
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in  such  a  way  as  to  develop  a  series  which  is  valid  both  at  the  turning 
point  and  away  from  it.  The  process  of  "piecing  together"  solutions  is 
thereby  avoided  and  the  proof  of  asymptotic  validity  becomes  relatively 
direct, 

I 

Now  it  will  not  be  surprising  to  find  that  the  nature  of  the 

present  problem  at  the  turning  point  depends  upon  whether  ^  is  in  Class 

2 

II  or  Class  III  and  in  the  latter  case  on  the  behavior  of  c^ri  at  x*. 

Thus  several  types  of  problems  involve  us  here  and  it  will  be  convenient 
to  employ  both  approaches  to  the  turning-point  problem.  Most  of  the 
problems  encountered  here  can,  in  fact,  be  solved  by  the  first  approach, 
without  the  necessity  of  attacking  the  connection  problem  at  all,  while 
the  Langer  method  seems  more  direct  for  the  purpose  of  completing  the 
analysis.  The  remainder  of  this  section  is  concerned  solely  with  the 
first  approach  and  requires  no  additional  analysis  beyond  the  application 
of  the  condition,  Eq,  (4.l). 

Class  II  Values  of  ^ 

For  the  sake  of  simplicity,  let  us  suppose  that  a  single  point  at 
which  u  =  ^  exists;  it  will  be  evident  that  the  procedure  described  here 
applies  equally  well  irrespective  of  the  nimiber  of  such  points.  Accord¬ 
ing  to  Sec.  3j  there  exist  solutions  0^,  ...  defined  on  <  x) 

and  a  second  set,  0^,  6^,  ...  defined  on  <  x  <  x*) .  Rather  than 

deal  with  these,  however,'  it  is  simpler  to  redefine  the  0^  to  be  the  set 
of  functions  obtained  from  the  set  valid  on  by  analytic  continuation 
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throughout  a  region  K  of  the  complex  z-plane.  Now  the  matrix  T(z, 
the  analytic  continuation  of  to  complex  is  differentiable  in 

some  strip  centered  on  the  x  axis^  punctured  in  the  neighborhood  of  x*^ 
as  showi  in  Fig,  1^  and  in  fact  is  single-valued  there.  Since  9^, 

Gg,  ...  involve  only  functions  which  are  regular  throughout  K,  and 
integrals  thereof,  it  follows  that  as  long  as  5C  is  simply  connected,  as 
in  Fig.  1,  0^,  Gg,  ...  are  regular  and  single -valued. 


Im(z) 


SHOCK 


Im(hj|)=CONSTANT 


1 

Uj 

J\ 

X 

,  >  CUT  C 

Re(h.,)=CONSTANT 

Fig,  1  The  region  K  of  the  complex  plane  through  which  the  series 
solutions  0^  can  be  defined  by  analytic  continuation  for  Class  II  values 
of  The  turning  point  x*  and  its  neighborhood  are  excluded.  By  virtue 
of  the  cut  C,  K  is  simply  connected.  The  light  curves  are  lines  of  con¬ 
stant  Re(hj2^)  and  the  paths  P.  of  the  type  illustrated  (heavy  curve)  have 
non-increasing  Re(h^P . 


4o 


We  can  now  investigate  condition  (4.l)  using  Eq.  (4.5).  Since 

and  are  real  and  positive  on  the  real  axis,  h^^  and  h^^  are 

monotonically  increasing  with  x.  In  region  K,  then,  the  lines  of 

constant  Re(h^.^)  are  perpendicular  to  the  x  axis,  as  shown  in  Fig.  1. 

Therefore  paths  P.,  originating  on  the  right,  can  reach  all  points  of  3C, 
0 

except  the  shaded  area,  simply  hy  proceeding  along  the  x  axis  and 
alternately  along  lines  of  constant  Re(h^.^)  and  Im(h^^)  with  decreasing 
Re(h.,),  as  shown  in  the  figure.  Since  this  argument  applies  for  both 

J-L  . 

j  =  2  and  5j  it  follows  that  9^  is  admissible  throughout  K,  ex;cept  in 
the  neighborhood  of  the  cut.  Thus  we  conclude 


0  ~  0^  for  all  X  ^  X*. 


(^.7) 


Class  III  Values  of  ^ 


Turning  points  where  c-  t]  increases 


In  this  case^  we  begin  by  considering  the  case  of  a  turning  point 
2 

through  which  c^t]  increases,  although  the  case  in  which  it  decreases  is 

related  by  a  simple  sign  transformation  of  the  independent  variable, 

2 

These  cases  in  >;hlch  c  ti  is  not  an  extremum  at  x*  can  be  reduced  to  an 

o' 

application  of  a  minor  generalization  of  the  "singular  turning-point 
15 

problem"  for  second-order  differential  equations,  the  generalization 
being  required  by  the  increased  number  of  equations,  i,e,,  n  +  4 
(rather  than  2)  equations. 

As  for  Class  II  values  of  we  begin  by  generalizing  the  intervals 
of  definition  of  the  series  solutions  to  the  complex  z-plane.  Now, 
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hoVever,  the  point  x*  is  a  branch  point  of  the  z-plane  analytic  continua,- 
tion  of  the  square  root  s(x,  C)j  whence  x*  is  also  a  branch  point 
of  and  Analytic  continuation  of  6^  and  will  yield  results 

which  depend  upon  the  path  of  continuation  and  some  care  must  be  exercised 
with  regard  to  notation. 

Begin,  then,  with  the  functions  0^,  6^,  6^,  ...  ,  defined  in  Sec.  5, 
on  the  interval  and  compute  M-22^(z)  and  in  the  neighborhood  of  x*, 

1^21(2)  =  [2/cd2/T)u]^^(z  -  x*)2[l  +  0(z  -  X*)] 


M-jl(z)  =  +  0(z  -  X*)] 


(4.8) 


where 

and  where  we  have  ta.ken  ^  to  lie  on  the  upper  imaginary  axis.  We  need 
not  consider  the  conjugate  points  at  all,  since  roots  of  V(t,  e)  occur 
in  conjugate  pairs,  as  is  proved  in  Appendix  B.  It  follows  from  Eq.  (4.8) 
that  for  z  «  X*, 


h2^(z)  -  h*^  =  [4Kd2/3iiuJ^(z  -  x*)5/^  +  ... 
hj^(z)  -  -  X*)  +  ... 


htj  =  h.  .(X.) 


The  loci, 

EeChjj^Cz)  -  hj^]  =  0 


(t.9) 
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Re[hj^(z)  -  h|^] 


0 


(4.10) 


are  referred  to  as  transition  lines,  whose  significance  will  be  evident 
in  what  follows.  Near  x*,  we  see  that  these  are  radial  lines  emanating 
from  X*,  viz. , 

cp  =  arg(z  -  X*) 

=  ±  rt/5j  ±  rt,  ±  5ir/5j  ••• 

^51  =  0,  ±  rt 

where  the  subscripts  classify  these  values  of  cp  in  accordance  with  the 

h„,  or  the  h,^  condition.  Note  that  the  line  cp  =  ±  n  and  the 
21  51  51 

cpg^  =  ±  Tt,  ±  5jt,  ...  lines  are  exact  transition  lines,  while  the  others 
hold  only  in  the  limit  z  -*  x*. 

Now  the  solution  9^,  which  is  defined  and  whose  validity  has  been 
established  for  x  on  presumably  will  retain  its  validity  in  the 
z-plane  near  To  examine  to  what  extent  this  is  true,  we  refer  to 

Fig.  2  which  shows  the  transition  lines  for  cp^^  and  the  sign  of 
Re(h2^  -  h^^)  for  the  neighborhood  of  x*  and  -  jt  <  cp  <  jt.  It  is 
observed  that  Re(h2^)  diminishes  on  paths  originating  at  the  right  and 
these  paths  can  be  extended  outside  the  sector  -  rt/3  <  cp  <  Jt/3>  to 
include  the  entire  plane  with  the  exception  of  the  axis,  cp  =  ±  jt,  by 
employing  paths  of  constant  Re(h2^)  as  shown  in  the  figure.  The  axis 
is  excluded  by  virtue  of  the  fact  that  0^  is  singular  at  s  =  0  and  hence 
no  path  can  be  extended  through  the  turning  point. 

To  proceed  with  the  description  of  the  asymptotic  solutions  of  the 
differential  equations,  we  now  introduce  some  additional  notation.  ¥e 
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Fig,  2  Paths  in  the  complex  z-plane  on  which  Re(h2^)  diminishes  a,nd 
Pj  on  which  Re(hj^)  diminishes  for  asymptotic  solution  9^,  defined  by 
analytic  continuation  from  the  real  axis  x  >  x*  to  points  with  |cp  |  <  rt, 

[cp  =  arg(z  .  X*)],  Points  z^  and  z^  on  opposite  sides  of  the  slit  at 
cp  =  ±  jt  both  lie  on  such  paths,  although  ^2(^2^  shown. 
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define  "sectors"  in  the  (z  -  x*) -plane  by  means  of  the  transition 
lines;  near  x*  these  approach  the  sectors 


Sq  :  -  rt/5  <  cp  <  rt/3 
;  rt/5  <  cp  <  rt 
S^;-rt<cp<  -  n/3 
Sg  :  rt  <  cp-  <  5jt/3 


(4.32) 


If  we  denote  the  analsrtic  continuation  of  solution  B .  to  sector  S  .  by 

J  ^ 

,  then  our  discussion  in  the  last  paragraph  proves  that  0^^^ ,  9^^^ , 
form  a  solution,  which  we  denote  simply  i.e,,  by  its  central 

sector,  asymptotically  admissible,  with  respect  to  condition  (4.l)  for 
hg^,  in  the  plane  slit  at  cp  =  ±  rt.  ¥e  refer  to  the  latter  as  the  Stokes 
lines  for  since  the  continuation  of  across  them  has  not  been 

shown  to  he  valid  (and,  in  fact,  is  not  valid). 

Similar  considerations  can  be  applied  to  the  solution  Gg  for  we 
have  near  x*. 


h^(z)  -  h^  =  -  [l|-/cd2/3T|u]^(z  -  x*)^/^-^  ... 

"  ^*2  “  -  X*)  H-  ...  (4.13) 

yielding  precisely  the  same  set  of  transition  lines  as  for  0^.  In  Fig.  3 

are  shown  the  transition  lines  for  h^g  on  -  rt/3  cp  <  5^/3  ^Jcid  the  sign 

of  Re(h  _  -  h*  )  is  shown  in  the  sectors  S_,  S  ,  and  Sg.  Paths  of 
12  12  u  X  ^ 
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Re(h.„-h.I)>0 


Re(h.«-h*  )<0 


9  STOKES  LINE 


Re(h,2-h,2)<0 


A''v' 


9  STOKES  LINE 


“^rt 

0  STOKES  LINE 


</)=5Tr/3V^<^=-Tr/3 


Fig.  3  Paths  P  in  the  complex  z-plane  on  which  Re(h_„)  diminishes  and 
P^  on  which  Re(h^2)  diminishes  for  asymptotic  solution  6'  defined  by 
6inal3rtic  continuation  of  02  from  the  real  axis  x  >  x*  to  points  with 
-  rt/3  <  cp  <  5rt/5,  Points  and  Zg  on  opposite  sides  of  the  slit  both 
lie  on  P^  paths  (that  for  Zg  is  not  shown)  but  two  types  of  P^  paths  are 
required.  The  dashed  P^  paths  cannot  enter  sector  Sg,  while  the  solid 
P,  paths  cannot  enter  the  cp  ^  0  region. 


decreasing  Re(h^)  then  originate  in  sector  and  cover  the  plane  slit 
at  9  =  -  rt/5,  5rt/3.  The  solutions  9^^^,  denoted  simply 

are,  therefore,  admissible  with  respect  to  condition  (4,l)  for  h^,  and 
have  Stokes  lines  at  cp  =  -  Jt/3j  5^/3. 

These  considerations  can  be  applied  to  the  continuation  of  the  0^ 
and  series  into  other  sectors  to  establish  a  collection  of  solutions 
which  are  admissible  on  certain  slit-planes,  as  seen  in  Table  I,  -where 
again  condition  (4»l)*has  been  applied  only  bet-ween  h^  and  h^.  Although 
this  list  could  be  extended  indefinitely,  it  is  readily  seen  that,  at 
each  point  z,  only  t-wo  distinct  solutions  exist.  Thus  and  0^^ 

differ  by  a  constant  factor,  so  that  Table  I  contains  but  three  distinct 
solutions,  say  G'’  %  B'  and  B'-  ,  although  solutions  B'  and  B' 

■will  also  prove  useful o 

TABIiE  I 

System  of  asymptotic  solutions  of  differential  Eq_.  (2,7b),  admissible 

with  respect  to  condition  (4,l)  for  (i,  j)  =  1>  2,  near  a  turning  point 

2  .  .  . 

at  which  c  ti  is  increasing, 
o 

Solution  _ Definition  Stokes  lines  a,rg(z  «  x^) 


^(0) 

#0) 

1  ’ 

^(l) 

^1  » 

3(-i) 

±  TC 

^(1) 

I— 1 

^(2) 

®2  » 

-  «/5, 

^(-1) 

«2  » 

e(-2) 

«2  > 

e(°) 

-  5rt/3, 

rt/5 

^(2) 

^(5) 

®1  » 

7rt/5 

^-2) 

^(-2) 
©1  > 

^-5) 
«1  > 

-  7n/3, 

-  rt/3 

4? 


It  is  of  interest  to  note  that  in  the  central  sector  of  the  plane 
slit  by  the  Stokes  lines,  the  solution  is  recessive  but  is  dominant 
in  its  outer  sectors.  For  example  in  S^,  f(0) 

is  recessive,  since 

hgi  -  h*^  has  positive  real  part  and  grows  in  magnitude  relative  to 

o)  -♦f  l) 

6^  'as  €  increases.  In  S^,  however,  Sg  'is  recessive.  In  fact,  it  is 

Just  this  property  of  the  central  sector  which  is  used  in  establishing 

the  condition  (4,l). 

Turning  now  to  the  validity  of  this  condition  for  J  =  3}  we 
employ  the  second  of  Eq,  (4,9)  and  (4.15)  from  which  it  is  seen  that 
Pj  paths  must  proceed  in  the  upward  direction  in  the  z-plane.  For 
it  is  seen  from  Fig,  2  that  such  paths  exist  for  all  points  in  the  region 
previously  found  to  be  admissible  with  respect  to  the  hg^  condition. 
Similarly  for  6^“^^  which  has  Stokes  lines  at  9  =  rt/3,  -  5rt/3j  the  slit- 
plane  is  entirely  covered  by  pa,ths.  For  however,  the  Stokes 

lines  extend  downward,  as  in  Fig.  5>  so  that  for  points  with  Im(z)  ^  0, 
the  Pj  paths  must  originate  on  the  same  side  of  the  Stokes  line.  Thus 
the  hj  condition  splits  this  solution  into  two  portions; 

valid  in  0  <  cp  <  53t/3 

9^^^  valid  in  -  rt/3  <  cp  <  rt.  ♦  (4,l4) 

Since  9'  '  differs  from  9^  '  by  a  factor,  a  similar  splitting  is  imposed 

on  it. 
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2 

Taming  points  where  c^t]  decreases 
2 

For  turning  points  at  which  c^tj  decreases,  the  solutions  and  their 
Stokes  lines  are  somewhat  altered.  It  is  straightforward  to  apply  the 
above  methods  to ‘obtain  the  corresponding  solutions.  We  begin  with 
0  ,  Op,  ...  ,  not  on  ,  but  to  the  left  of  x*  and  proceed  by  analytic 
continuation  as  before.  Thus,  the  first  of  Eq.  (4.9)  becomes 

hp^Cz)  -  h|'^  =  [V/c(-d)2/5Tiu]^(x*  -  +  ...  (4.15) 

Defining  sectors  S^,  S^,  ...  by  the  transition  lines  which  near  x*  are 
given  by 

5p  =  arg(x*  -  z) 

$21  =  ^  ±  ••• 

$51  =  0,  ±  rt 

as  indicated  in  Fig.  4,  we  can  readily  obtain  the  system  of  solutions  in 
Table  II.  The  regions  of  validity  refer,  as  in  Table  I,  only  to  the 
condition  (4,l)  applied  between  6^  and  Gg,  Note  that  Og  is  now  recessive 
in  Sq. 

Imposing  the  condition  (4,l)  for  j  =  5  leads  again  to  paths 

7?f-l)  ”*(2) 

directed  upward  in  the  z-plane^  whence  0^  and  6^  are  seen^  from 
Fig,  4,  to  be  split  into  two  parts,  for  example; 

valid  in  -  5i^/5  <  9  0 

valid  in  -  rt  <  9  rt/5»  (4,l6) 
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Fig.  4  Definition  of  sectors  S^,  S^,  and  S  ^  for  a  turning  point  x* 
through  which  c^t]  decreases.  The  sign  of  ReCh^^  -  h2^(x*)]  obtained 
from  the  analytic  continuation  of  and  0^  from  the  real  axis  x  <  x* 
is  indicated.  paths  on  which  condition  (4,1)  is  valid  proceed 


upward  in  the  z-plane,  as  indicated. 


TABLE  II 


System  of  asymptotic  solutions  of  differential  Eq,  (2,2),  admissible 

with  respect  to  condition  (4,l)  for  (i,  j)  =  1,  2,  near  a  turning  point 
2 

X*  at  which  c^t]  is  decreasing. 


Solution 

Definition 

5(0) 

■5(0) 

eh) 

^2  » 

I— 1 

5(0) 

o[-^K 

5(0) 

5(2) 

5(2) 

“2  > 

qO) 

®2  » 

5(-2) 

5(^) 

«2  » 

0(-5) 
^2  » 

-  «/5,  5«/5 

-  5rt/3,  Tt/5 

rt/3,  7rt/3 

-  7rt/3,  -  rt/5 


Continuation  problem 

The  continuation  of  the  asymptotic  expression  for  0  beyond  a 

turning  point  divides  into  t™  problems,  then,  depending  upon  ‘whether 
2 

c^Tj  is  increasing  or  decreasing  at  x*j  and  a  complete  determination  of 

0  for  all  X  will  depend  upon  the  configuration  of  turning  points  along  the 

•  ^ 

X  axis.  It  is  convenient  to  carry  out  the  determination  of  0  under  the 
2 

assumption  that  c^t|  has  some  simple  form,  although  more  complicated  forms 
can  be  treated  with  little  additional  difficulty.  Three  profiles  will  be 
considered,  as  shown  in  Fig.  5j  each  occvirs  for  certain  values  of  the 
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parameters  for  the  idealized,  one-reaction  detonations  whose  stability 

k 

has  been  studied  numerically.  Profile  I  designates  those  steady  detona- 

2  2 
tions  in  which  c^ri  increases  monotonicallyj  profile  D,  those  in  which  c^ti 

2 

decreases  monotonic ally;  and  profile  M  designates  those  in  which  c^ti  has 
a  single  maximum  as  a.  function  of  x. 


2 

Fig.  5  Profiles  of  c  t)  in  the  steady  flow  which  are  considered.  Profile 

°  2 
M  includes  both  cases  with  regard  to  the  relative  values  of  c^t]  at  x  =  0 

and  X  =  “  , 
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Profile  I,  For  detonations  ha,ving  profile  I  behavior,  we  first  note 
that  (for  ^  in  Class  III)  to  the  right  of  x*,  i.e.,  on  interval  we 
liave  e  Eq.  (4.6),  and  c^ti  increases  at  x*.  The  asymptotic  expres- 

sions  of  Tatle  a.s  modified  by  Ecj.  (4*l4)^  a.re  directly  applicable# 

Thus  6^  is  contained  in  9^^^ which  brea.ks  down  at  cp  =  ±  jt .  To  solve  the 
continuation  problem^  we  write  0  as  a  linear  combination  of  two  independ¬ 
ent  solutions j  valid  on  cp  =  it  [see  Table  I  and  Ecj#  (4.l4)3 

and  on  cp  =  -  it 

0  =  020^“^^  +  (4.18) 

Contributions  from  the  remaining  n  +  2  solutions  need  not  be  included, 
for  reasons  which  will  be  given  later.  Now  in  6^  ^  and  Eq.  (4.17) 
are  valid.  Since  0^^^^  and  are  identical,  viz.  in  S^  and 

|(11)  is  recessive  there,  it  follows  that  d^  =  1.  Similar  considerations 
in  S  ^  prove  that  d^  =  1. 

The  remaining  coefficients  are  determined  by  the  condition  that  0 
is  single-valued.  Equating  Eq.  (4.1?)  and  (4.l8)  for  x  <  x*  and 
utilizing  the  fact  that  and  differ  by  a  constant  factor,  as  do 

0^^^^  and  'w®  s®®  that 

^  ^  e(-2i)  +  ^(2)  ^  e[-i)  +  (4.19) 
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That  is,  is  obtained  by  adding  the  analytic  continuations  of  9 ^ via 

both  the  upper  and  the  lower  half -planes. 

The  absence  of  the  n  +  2  other  solutions  from  Eq,  (4.17)  and  (4,l8) 
is  justified  by  the  consideration  of  6  in  S2  and  respectively.  Since, 

for  these  z,  h,^  and  h_o  have  larger  real  parts  than  at  x*,  any  contri- 

’  51  5^ 

button  of  a.  p^-type  solution  would  be  dominant  over  the  and  p^-type 
solutions.  Since  0  itself  is  known  asymptotically  in  S_^,  viz., 

,  it  follows  by  single-valuedness  that  p^-type  contributions  are 
not  permitted. 

•It  can  be  readily  seen  that  0^^^  contains  for  the  square  root  s, 

for  X  <  X*,  s  =  i  |s  I  and  hence  is  a  p^-type  solution,  while  0^”^^  has 

.  s  =  -  i  |s  I  and  hence  is  a  p^ -type  solution.  To  express  0  in  most 

convenient  form,  we  now  specify  0^  and  completely  (to  highest  order) 

by  specifying  the  lower  limit  of  integration  in  Eq.  (5.6)  and  (5.10)  to 

be  X  =0.  The  k.  integrals  then  exist  for  all  x  5/  x*  and  continue  to 
o  1 

have  derivatives  E. provided  the  integration  path  for  x  >  x*  avoids  the 

IX 

turning  point  by  a  small  excursion  into  the  upper  half  z-plane.  With 
this  understanding,  then. 


(4.20) 

♦ 

e  ~  +  oSg ,  for  x  <  X* 

with  the  coefficient  C(  found  to  be 


54 


-o)  =  exp  p(e,  -o) 


P(e,  v>) 


[eiJ-i  +  E^^ldz 


■where  C  is  the  contour  enclosing  x*  shown  in  Fig.  6. 


(4.21) 


Im  (z) 


Fig.  6  Contour  C  in  the  complex  z-plane  used  for  the  definition  of 

0)5  Eq.  (4,21)5  required  for  the  continuation  problem  at  a 

2 

turning  point  x*  at  which  c^t]  is  increasing  with  distance  from  the  shock. 
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2 

Profile  D.  For  Class  III  values  of  decreases  at  the  turning 

point.  The  a,symptotic  expressions  of  Table  II,  modified  by  Eq.  (4,l6) 

are  now  applicable.  For  x  >  x*,  we  had  by  Eq.  (4.6),  9  ~  e]_.  Now  the 

-»(l) 

entry  in  Table  II  -which  is  valid  for  x  >  x*  and  of  ki^-type  is  0^  ' ,  as 
is  readily  verified  by  consideration  of  in  the  neighborhood  of  x*. 

Now  S^^^is  contained  in  0^^^  which  is  valid  for  x  <  x*  as  well,  i.e., 
for  9=0.  Thus  for  profile.!). 


0^^^  for  X  <  X* 

l) 

0^  '  for  X  > 


If  we  now  utilize  definitions  of  6^^  and  Og  analogous  to  those  for  profile 
I,  viz.,  x^  =  0  in  Eq.  (5.6)  and  (5.10),  with  the  integration  paths  for 
h^,  hg,  k^,  and  a.voiding  x*  by  a  lower  half  z-plane  excursion,  this 
can  be  written  simply  as 


0  ~  0^  for  all  X  /  X*, 


(4.22) 


Profile  M.  To  treat  profile  M,  we  consider  the  value  of 
i 

^  =  i  (c  q^)  separately,  for  the  turning  point  in  that  case  apparently 

O  I11Q.X 

differs  from  those  we  have  already  considered.  We  return  to  its  study 

in  Sec.  6.  For  all  other  Class  III  values  of  we  simply  combine  the 

results  for  profile  I  and  profile  D.  We  denote  by  x*  the  turning  point, 

2 

if  any,  where  c^q  is  decreasing  and  simply  use  x*  to  denote  the  point, 

2 

if  any,  at  which  c^q  is  increasing.  For  those  ^  for  which  only  one 
turning  point  exists,  the  considerations  of  profile  I  or  D  apply.  For 
those  ^  for  which  two  t-urning  points  exist,  the  outer  one  x*  is  of  the 
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(4.25) 


0,  and 
and 

therefore 

(4.24) 

not  X*. 

The  asymptotic  evaluation  of  V(t,  e)  requires,  according  to 
Eq,  (2,7a),  the  determination  of  the  function  L(e,  u)  which  to  hipest 
order  in  £  is 

L(e,  x>)  s=  +  ih2  -  ©(O,  Cj  +  ih^) .  (5.l) 

The  evaluation  of  and  b^  will  evidently  be  the  major  task  of  this 
section,  since  the  scalar  product  is  essentially  known  from  our  expres- 
sions  for  9. 

^  in  Class  I 

Here  we  distinguish  between  those  ^  for  which  h^(x,  0  is  pure 
imaginary  and  those  for  which  it  has  a  real  part.  The  former,  Re(0  =  0, 

i 

lni(0  ^  (c  Ti^)  ,  we  denote  as  Class  laj  and  the  latter,  Re(C)^0> 
o  inQiX 

lm(0  ^®(^)  ^  ^  ^  Class  II,  we  designate  as 


profile  D  type.  The  continuation  problem  there  yields 
^  ~  "9  for  X  >  X*,  X  X* 

where  and  9-^®  defined  as  for  profiles  I  and  D,  i.e.,  x^  = 
the  h.,  k^  integration  paths  excise  x*  via  the  upper  half -plane 
X*  via  the  lower  half-plane.  The  continuation  problem  at  x*  is 
typical  of  profile  I,  and  we  obtain 

with  Oi  given  in  Eq.  (4,21)  and  with  contour  C  encircling  x*  but 

5.  ASYMPTOTIC  EVALUATION  OF  V(t,  e) 
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Class  Ib.  Now  the  integrals  h^  and  hg  "the  asymptotic  form 


b 


1 


o)  exp  Ceh^(x,  0]dx 


A(x,  -o)  =  exp  [k^(x,  u)]  [£(x,  0  *  e^l  *  A"^  *  (5.2) 

with  h^  and  A  regular  functions  of  x.  For  Class  la,  h^  is  pure 
imaginary  and 

^x 

\\  1  =J^  Odx' 

increases  monotonically  with  x.  Thus,  we  can  transform  to  z  =  |h^  1  as 
integration  variable  to  obtain  a  Fourier-type  integral,  viz,. 


/CO 

A(x,  C,  b)[ip  (x,  OT^dz  . 

Since  g^  is  proportional  to  the  chemical  rate  (see  Ref,  5)  which  vanishes 

exponentially  with  x,  as  x  -*  ®,  and  since  z  is  linear  in  x  for  large  x, 

it  follows  that  A  and  all  its  derivatives  with  respect  to  z  vanish  at 
A  l4 

z  =  “,  Thus  b^  can  be  evaluated  asymptotically  in  €  to  satisfy  the 
order  relation 

b^  =  0(l/e)  .  (5.5a) 

Identical  considerations  show  also  that 

bg  =  0(l/e)  .  ’  (5.5b) 

For  Class  Ib,  h^  has  a  real  part  which  is  negative  and  decreases 
linearly  with  x  as  x  -»  »,  Thus,  Eq,  (5.2)  yields  a  Laplace-type 
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integral, 

00 

kl  |a(x,  ?,  «)  I  [-Ee(n,)r^da 

-L  ^0 

z(x,  C)  =  -  Re[h^(x,  ?)]  . 

Since  |a  |  [-  Re(ii^)]"^  is  bounded  (say,  by  M)  ,  it  follows  that 

|bj  ^  yij  e“^.^dz  =  0(l/€)  .  (5.4) 

0 

The  same  order  relation  is  satisfied,  of  course,  by  b^. 

^  in  Class  II 

For  those  (real)  ^  in  Class  II,  ~  0^  remains  valid  except  in 
the  neighborhood  of  the  turning  points  x^.  To  evaluate  b^  and  ^2  for 
this  case,  the  path  of  integration  can  be  deformed  into  the  uppei*  half 
z-plane  so  as  to  excise  the  turning  points,  without  altering  the  values 
of  b^  and  b^.  Along  such  a  path  K,  ^  ~  6^  is  valid,  and  we  have  to 
highest  order 

b^  exp  [eh^(x,  A(x,  ^,'D)dx, 

as  in  Eq.  (5.2);  K  extends  to  x  =  ®,  of  course.  It  is  seen  from  the 
expression  for  E^^  in  Appendix  A  that  the  integrand  is  regular  at  the 

+ 

The  symbol  z  is  used  repeatedly  within  this  section  as  a  variable  of 
integration,  but  its  definition  varies  with  each  usage. 
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turning  points  (although  correction  terms  in  0^  might  be  singular  there), 
and  the  path  K  can  be  replaced  by  the  real  axis  once  more.  Just  as  for 
C  in  Class  Ib,  then, 

b^  =  0(l/€) 

b2  =  0(l/e).  (5.5) 

^  in  Class  III 

The  discussion  in  Class  III  becomes  considerably  more  complicated 
by  virtue  of  the  singalarities  in  and  6^  at  the  turning  point.  As 
a  preliminary  to  our  discussion,  it  is  useful  to  first  write  these 
asyn5)totic  expressions  0^  and  0^  so  that  the  infinity  associated  with  x* 
is  evident.  From  the  expressions  for  and  E^2  (which  determine 
and  kg)  in  Appendix  A,  it  is  seen  that,  for  turning  points  at  which  c^t] 
is  not  an  extremum,  k^^  and  kg  are  infinite  at  x*  only  in  the  contribution 
from  the  -  -gd  An  s/dx  term,  Eq.  (A.16);  the  l/s  terms  lead  to  finite 
contributions.  Thus  we  write 

k^(x,  C,  ^)  =  -  i  I'j  +  l#(x,  C,  ^) 

/X 

^  Xn  s/dx^]dx^  (5*6) 

0 


andj  using  subscript  +  to  denote  evaluation  behind  the  shock,  we  write 


Before  considering  "the  specific  evaluation  of  tlie  b  integrals  j  it 
is  profitable  to  make  some  general  observations  concerning  them.  The 
contribution  to  the  b's  from  the  neighborhood  of  the  turning  point  wild, 
require  special  attention,  for  the  asymptotic  expressions  of  the  last 
section  are  singular  at  a  turning  point  and  are  not  valid  in  its  neighbor¬ 
hood.  The  method  of  evaluation  is  that  introduced  for  Class  II  values  of 
the  path  of  integration  is  deformed  to  bypass  the  turning  points  with¬ 
out  affecting  the  values  of  the  integral.  A  certain  6 -neighborhood  can 
be  avoided  by  writing,  in  schematic  notation, 

/•x*-6  rO  ■im  -im 

b  =  /  l(x)dx  +  i6  /  l(x^  +  6e  ^)e  ^dcp 

*'0 


(5.7) 


The  integrand  l(x)  is  evaluated  then  from  the  asymptotic  formulas  arid  is, 
of  course,  singular  only  at  x*.  Thus,  the  result  for  b  will  be  inde¬ 
pendent  of  6  and  we  are  free  to  enpiloy  the  6-^0  limit  without  regard  to 
the  fact  that  the  asymptotic  expression  for  6  loses  its  validity  near  x*. 


Profile  I 


To  evaluate  b^  for  profile  I,  we  consider  four  contributions.  First 
from  the  (O,  x^)  interval,  we  have,  in  view  of  Eq.  (4,20),  contributions 


(5.8) 


6l 


Wow  and  i-i^  are  pure  imaginary  on  this  interval  and  the  branch  points 
of  and  ^2  ^  factored  by  a  change  of  integration  variable. 
Thus,  we  introduce  a  real  variable 


z  =  -  xs 

so  that  Eq.  (5»8)  becomes 


^11  ~ 


=  /  -"axp  [ie(Pj^(z)] 

A 

b^  =  y*  "^exp  [iecp2(z)]  z2A2(z)dz 


fz  2zp 

i^(z)  =  I®  /  ^  jt  ”3! — 

2(_  is^)^ 


A.(z)  = 


(5.9) 


2 

where  d  is  the  derivative  of  c^q  introduced  in  Eq,  (4.8),  and  where  the 
dependence  on  ^  and  o  has  been  suppressed.  Clearly  cp^(z)  and  A^(z)  are 

analytic  for  all  0  iS  z  s  so  that  we  can  apply  the  method  of  stationary 
15 

phase  to  obtain 

=  0(e-5/‘‘) 

b^  =  0(e-5/^).  '  (5.10) 

From  the  path  around  the  turning  point  arises  the  contribution 

h5  • 

'Jt 
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Since  the  integrand  is  of  order  6“^  for  small  6,  it  follows  that  in  the 
limit  6  -►  0, 

^3  =  0  (5.11) 


Finally  the  contribution  beyond  the  turning  point  is 


and  i-L^  is  complex  with  negative  real  part.  A  suitable  asymptotic  bound 
on  b^j^  is  obtained  for  its  absolute  value, 


A"^  ♦  g^  exp  (k*)  I  dx 


z  =  (K:s/T)u)dx^, 


Now  it  can  be  seen  that  z  has  a  power  series  expansion  in  s, 
beginning  with  the  cubic  term.  Thus  conversion  to  the  integration 
variable  z  yields 

l^lJ  e"®^z"2A^(z)dz 

Jl 

^  •  y  •  •  St  • 

Since  g^  vanishes  exponentially  in  x  and  1#  has,  for  Re('o)  ^  0, 
negative  real  part  and  is  linear  in  x  [Eq,  (5»8)],  it  is  seen  that  is 
bounded,  whence 
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(5.35) 


KJ  ^  m/*  e-^V^dz  =  0(€‘2). 

To  establish  the  order  of  b^,  then,  we  yet  require  the  order  of  the 

coefficient  Of  of  b^.  But  Of  depends  on  e  only  through  the  integral  of 

p  (i.e.,  h  )  which  is  pure  imaginary  on  the  x  <  x*  portion  of  contour  C 
3L  1 

of  Eq..  (4.21).  Since  is  finite  at  xJ^  itself,  the  portion  of  C  around 
X*  yields  zero  in  the  limit  of  vanishing  radius,  whence  we  obtain 

«=0(e°)  (5*15) 

and  we  combine  Eq.  (5 .10)  >  (5  .H)  j  (5 .12)  j  and  (5 .15)  to  yield 

b^=0(e“2).  (5.14) 

Identical  considerations  apply  to  b^  as  well. 

Profiles  D  and  M 

The  derivation  of  order  relations  in  both  of  these  cases  is  similar 
to  that  for  profile  I  and  in  both  cases  the  relation  (5.l4)  is  recovered. 

Expressions  for  L 

It  follows,  then,  that  the  b’s  can  be  ignored  in  Eq.  (5.l)  and  only 
the  scalar  product  remains.  For  ^  in  Class  I  and  Class  II,  the  only 
contribution  to  ?(0,  n,  e)  is  that  from  which  yields 
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L(e,  i,  ~  L^(0 


u  (1-v  )  +  W  \1 

L  fO  =  -  r  I  ■  I  7.U.  . . —  +  T]  ( 1  -  '  - ■  —  I 


(5.15) 


i  =  2  -  (1  -  n)(i  -  %)v_Ps/t 


The  new  symbols  in  Eq.  (5.15)  are  Tj  the  temperature,  and  Pg,  the 
entropy  derivative  of  pressure  p  at  constant  volume  v  and  composition 
Here  we  have  used  subscripts  +  and  -  to  denote  evaluation  in  back 
of  and  in  front  of  the  shock,  respectively.  This  is  simply  the  value 
of  V  for  a  step  shock.^ 

For  ^  in  Class  III,  the  value  of  9  at  the  shock  depends  upon 
whether  we  are  considering  profile  D  or  profiles  I  and  M,  In  the 
former  case,  we  still  have  0  ~  0^  so  that 

L  ~  Ij^(0  profile  DJ  (5.l6) 

2 

but  for  those  cases  where  c^tj  increases  with  distance  at  the  shock,  we 
have 

L  ~  L^(0  +  ^{^3  ?>b)L2(0>  profiles  I  and  M  (5 .17) 


where  1^(0  differs  from  L^,  Eq.  (5. 15),  only  in  the  appearance  of  -  s_^ 
instead  of  s^,  i.e., 

u  (1-v  )  P  ^Uk  s  )  /  .2  v-j 
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Detemination  of  Of 


The  coefficient  a  of  the  non-shock  (L^)  contribution  to  L  for 
profiles  I  and  M  is  defined  in  Eq.  (4.2l).  Letting  the  radius  of  the 
circulax  part  of  C  approach  zero,  one  finds,  using  ^ from  Appendix  A, 
the  logarithm  of  a  to  be 


P(6,  C,  =  Xn  a(e,  C, 

P(e,  C,  ■o)  =  rti/2  -  ieP^(^j.)  +  ^2^^^  “ 


x*(C^) 


2K  s  dx 


Tlu 


,  .  -c  r  1 

2  -  Tp 


r  1  dT|  a*r 

[_1-T1  dx  "  T  u 


2u?»? 

(i-n)(u2+t2)  ■  11  J 


dx 


^  / 

*D 


2/< 


nls  I  u 


dx 


(5.18) 


The  subscript  v  on  the  rate  r  denotes  the  partial  derivative  at  constant 
entropy  and  ccmposition.  We  note  that  Pg,  and  are  all  real  and 
that  p^  and  P^  are  positive  but  Pg  depends  on  the  equation  of  state  and 
kinetics  of  the  particular  system.  For  profile  I,  however,  it  is  seen 
that  Pg  becomes  negative  for  sufficiently  near  the  value 
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aM  in  fact  approaches  -  oo,  by  virtue  of  the  fact,  noted  in  Sec,  5j  that 
only  the  a  •  r_^  terra  continues  to  contribute  to  the  integral  for  large  x, 

^4  -1 

and  a  •  r^  is  necessarily  positive  at  chemical  equilibrium. 

Stability 

At  the  outset,  it  is  expeditious  to  assume  that  the  equation  of 

state  of  the  unreacted  explosive  is  such  that  the  von  Neixmann  shock  is 

Itself  stable  (i.e.,  it  would  be  stable  if  the  reactions  behind  it  were 

somehow  suppressed) ,  for  we  thereby  concentrate  attention  on  instability 

which  arises  solely  from  the  reactions.  Thus,  we  assume  that  tias 

no  roots  in  the  right  half -plane,  and  in  fact  we  require  that  its  roots 

have  negative  real  part.  It  follows  that  the  asymptotic  expression  for 

V  is  non-zero  for  all  Class  I  and  Class  II  values  of  whence  any 

izns table  roots  are  to  be  found  with  ^  pure  imaginary  and  in  Class  III 

with  Re('D)  >  0,  .Moreover,  detonations  characterized  by  profile  D  are 

stable  for  large  e,  since  the  shock  expression  for  V  holds  for  Class  III 

values  of  ^  as  well.  Thus  we  now  confine  our  attention  to  profiles  I  and  M, 

Wow  the  behavior  of  has  been  given  in  the  discussion  of  shock 

stability,^  from  which  it  can  be  seen  that  our  assumption,  that  L^(C) 

only  has  roots  with  Re(^)  <  0,  implies  that  li^(0  is  positive  for 
2  2 

and  increases  monotonically  with  (see  Case  Ib  of  Ref.  6). 

Wow  cc,  Eq.  (4,21)  and  (5»l8),  is  periodic  in  €  and  has  magnitude 

|a  I  =  exp  [pg  -  Re('o)p^]  . 

If,  therefore,  for  ^  in  Class  III,  we  have  the  inequality 
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(5.19) 


it  follows  that  for  Re('D)  ^0,  L  =  +  CCI^  has  positive  real  part 

irrespective  of  €  and  Im('o) ,  Hence,  for  those  for  which  Eq.  (5  •19) 
holds,  V  has  no  unstable  roots. 

However,  if  for  some  we  have 

|l^(iy  I  exp  [^2(^1  (5.20) 

then  unstable  roots  occur  for  Re(u)  satisfying 

=  |l^  I  exp  [^2  -  Re(u)p^]  (5.2l) 

at  periodically  distributed  values  of  Im('o),  viz. 

ImCu)  =  [(2n  -  Dit  -  >  0, 

Im(i3)  =  [(2n  +  <  0,  (5.22) 

for  n  =  0,  ±  1,  ...  .  It  follows,  then,  that  every  point  on  the  vertical 

line  in  the  u-plane  satisfying  Eq.  (5 .21)  is  asymptotically  a  root  of  V 

for  periodically  distributed  values  of  e, 

e  =  [(2n  -  |-)it  -  Im(o)p^]p"^,  >  0, 

e  =  [(2n  +  |)it  -  Im(o)p^]p“^,  <  0,  (5.25) 

for  integer  values  of  n  sufficiently  large. 
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To  understand  the  behavior  of  these  roots,  we  consider  the 

asymptotic  expression  for  V  for  finite  e.  Let  us  suppose  that  L^/L^ 

2 

and  exp  (^2)  behave  as  in  Fig.  7j  which  is  appropriate  for  a  c^t)  profile 
I,  The  value  zero  at  the  ri^t  for  exp  (p^)  occurs  by  virtue  of  the 
approach  of  to  -  °=>  as  the  turning  point  recedes  to  x*  =  Between 
^il  ^i2’  unstable  roots  occur,  as  given  above,  while  at 

^il  ^i2  stable  roots  occur.  This  interval  is  seen  to 

correspond  to  a  series  of  unstable  €-intervals,  •which  for  Im('o)  =  0  are 


Fig.  7  Sketch  of  the  functions  and  exp  (p^)  for  profile  I 

(see  Fig,  5)  detonations  for  =  Im(T)/€  traversing  the  Class  III 
segment  of  the  ^-plane.  Instability  is  obtained  in  the  e  -*  «>  limit 
for  a  detonation  which  has  exp  (Pg)  greater  than  |l^/i^|  ,  as  in  the 
interval  . 
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for  n=l,  2,  ,,,  .  Asn  increases  these  intervals  overlap,  so  that  for 

e  above  some  value  all  wave  nvimbers  are  unstable  and  no  estimate  of  the 
range  of  unstable  e  values  is  obtained.  If  the  exp  (Pg)  ^1^^  curves 
do  not  cross,  then  no  unstable  roots  occur  for  large  e. 

Two  additional  points  should  be  noted.  First  we  observe  that 

2  2  2  2 
and  are  equal  at  Since  we  had  Of  =  0  at 

for  profile  I,  it  follows  that  the  asymptotic  expression  for  V  behaves 

continuously  as  traverses  the  boundaries  between  Class  I  and  Class  III 

for  profile  I  as  well  as  profile  D.  Because  of  this  continuity,  the 

2  2 

asymptotic  value  of  V  for  =  ,  a  point  excluded  from  all  our 

considerations,  is  taken  simply  V  =  eL^(0>  as  is  obtained  by  approach 

to  this  point  from  either  Class  I  or  Class  III. 

A  second  point  for  considerations  is  concerned  with  profile  M.  It 

2  2 

is  evident  that  V  behaves  continuously  with  both  at  =  (Vi  and 

2  2  2 
^.  =  (011)  •  At  the  maxirrrmn  value  of  c  ti,  however*  continuous  behavior 

can  be  lacking.  Although  ^oes  become  infinite  at  this  point,  as  can 

be  seen  from  Eq,  (5  #18),  its  sign  can  be  either  +  or  «  simply  on  the  basis 

of  the  kinetic  coefficients  r  •  If  0£  is  also  infinite,  whence 

2  2 

dominates  L^,  while  for  only  the  term  appears  in  the 

asymptotic  formula.  It  is  clear,  then,  that  the  complete  stability  be¬ 
havior  of  profile  M  requires  an  asymptotic  expansion  for  0,  valid  at  the 


maximum.  This  task  is  reserved  for  the  next  section 


6.  TUBKING- POINT  PROBLEM  FOR  AN  EXTREMUM  IN 


Althou^  the  method  of  Sec.  4  can  presumably  be  employed  to  solve 
the  continuation  problem  facing  us  here,  it  is  certainly  a  far  more 
difficult  problem.  In  particular,  the  discussion  of  Ref.  5  indicates 
the  need  to  solve  a  connection  problem  of  a  rather  difficult  sort. 
Instead,  we  turn  to  the  so-called  "related-equation  method"  of  Langer, 
with  the  hope  that  ranch  of  the  present  problem  has  already  been  solved 
in  its  essentials  by  virtue  of  the  known  solution  of  certain  turning- 
point  problems  for  second-order  equations.  These  second-order  equations, 
of  the  form 


d^u 

dx2 


^aPq^Gx)  +  a5q^(x)  +  R(x,  a!))J  u  =  0 


R 


=2 

i=0 


(6.1) 


(the  absence  of  a  term  in  du/dx  implies  no  loss  of  generality^),  have 
been  solved  asymptotically  in  m  for  the  case  of  a  simple  zero^  of  <3^(x) 


and  a  second-order  zero.^^  By  writing  Eq.  (6.l)  as  a  systan 


u^  =  u 


du. 


rs —  = 
dx  2 


du/dx  =  OD^  •  u 
0 


(6.2) 


Q  = 


-  (q^  +  +  R/' 


-^2 

CO 


1 

0 
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we  see  that  the  simple  zero  corresponds  to  the  problem  of  Sec,  4  in  so 

far  as  the  difference  in  the  eigenvalues  of  the  leading  matrix  vanishes 
1 

as  (x  -  x*)^,  while  the  case  of  a  second-order  zero  of  coincides, 

in  the  same  sense  with  the  problem  at  hand.  The  present  section  is 
based  upon  the  exploration  of  this  correspondence. 

The  procedure  which  will  be  followed  is  as  follows.  Non-singular 
transformations  of  our  equations  are  devised  which  yield  a  second-order 
system  of  the  form  of  Eq,  (6,2)  and  an  (n  +  2) -order  system,  with 
interactions  between  the  two  systems  first  appearing  with  terms  of  order 
in  the  coefficient  matrix,  for  arbitrary  m.  The  related  equations 
are  taken  to  be  independent  systems  of  order  2  and  order  n  +  2,  with 
coefficient  matrices  agreeing  with  the  transformed  system  also  to  order 
e"*^.  The  significance  of  the  related  equations  derives  from  the  fact 
that  their  solutions  are  known.  In  particular,  the  related  equation  of 
order  2  is  taken  to  be  the  related  equation  of  McKelvey,^^  while  the 
asymptotic  series  given  by  Friedrichs  form  the  related-equation  solutions 
for  the  (n  +  2) -order  equations.  The  proof  that  there  are  actual 
solutions  asymptotic  to  the  related-equation  solutions  has  essentially 
been  given  by  McKelvey  and  is  not  detailed  here. 

Related  Equation 

We  begin  once  more  with  Eq,  (2,7b) 

[6^(x)  +  ^(x)]  •  e  (6,3) 
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and  restrict  our  attention  henceforth  to  Class  III  values  of  ^  and 
suppose  that  the  x  interval  under  consideration  contains  a  single 
turning  point.  We  -wish  to  obtain  a  related  equation,  for  any  integer 


m  ^  1, 


dcp 


00 


(6.4) 


such  that  a  complete  system  of  n  +  4  solutions  of  Eq,  (6,4)  are  known. 
To  accomplish  this  task,  we  proceed  inversely  by  transforming 
Eq,  (6,5)  so  that  it  will  resemble  equations  with  known  solutions.  We 
look  for  a  similarity  transform 


e  =  So(^>  •  ^(0) 

^  axp  +  H2)ax'J 


(6.5) 


such  that,  to  hipest  order  in  the  coefficient  matrix,  the 
differential  equation  will  be  composed  of  a  second-order  equation  of 
the  form  Eq,  (6,2)  and  a  diagonal,  (n  +  2)-order  system.  In  the 
equations , 
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-3^  =  ]  •  *(0) 

•  [ib  -  i  (“i  *  “a’s.]  •  & 

IM'S’’"  [4'£o-'^V'^] 


(the  dependence  on  ^  and  n  being  suppressed  for  the  present)  we  observe 
that  ^  has,  in  view  of  Eq.  (5.5)  j  eigenvalues,  ^  M.^,  '2^21’ 

-  ^22^  whence  we  regiiire  ^  be  given  by 


ilo 


-Xq(x)  0 


0 


(6.7) 


Xq(x)  =  -  ^2^/4  =  -  (1  -  Ti)(^^  +  c^ti)/ti^u^ 


n(x)  =  “  ^^2l/2  = 

where  I.  denotes  the  unit  matrix  of  order  i.  It  is  important  to  note 
that  the  exponential  in  Eq,  (6.5)  and  the  matrix^  do  not  contain 
algebraic  branch  points  at  a  turning  point,  for  the  sqaare  root  s  does 
not  appear.  Thus,  ^(x)  can  be  found  as  a  regular  function  of  x,  and 
this  determination  is  given  in  Appendix  C,  In  addition  ^  is  a  non¬ 
singular  matrix  (i.e,,  1^|  ^  O)  for  all  x.  In  view  of  these  properties 
of  it  is  clear  that  H(x),  Eq,  (6.6),  is  a  regular  function  of  x. 
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To  proceed  further,  we  look  for  additional  similarity  transforms 
which  will  turn  the  lower-order  matrices  of  the  differential  equation 
into  block-diagonal  form,  with  the  upper  2-by-2  of  the  special  form 
displayed  by  Eq,  (6,2),  and  at  the  same  time  not  affect  the  leading 
matrices  which  are  already  in  this  form.  The  existence  of  such  matrices, 
of  the  form 

r  g,.  ,  (6-8) 


l6 

for  k  ^  1,  is  suggested  by  the  work  of  Turrittin,  We  let 


e  =  s 


U) 


(6.9) 


and  require  that 


dt 


jM  .  64(x,  €)  . 


X,  e)  =  €)  + 


c)  = 


(6.10) 


u 


(k) 


U^^^(x) e"^ 

i«-V  '  ' 


L  (x)  = 


0 


JL  (x) 
v'  ' 


0 

0 


,  V  ^  1  . 
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The  proof  that  such  exists  is  also  given  in  Appendix  C  and,  again, 
these  transforms  are  regular  functions  of  x.  Obviously  each  has  an 
inverse  for  e  sufficiently  large. 

For  any  given  integer  m,  the  transformations  in  Eq,  (6.9)  can  be 
terminated  to  yield 


?(x,  e)  =  e) 

00 

U(x,  e)  =  U^”^\x,  e)  =  2!^ 

^(x,  e)  =  e) 

^(x,  €)  =  e) 

[€^(x,  e)  +  e'^^VCx,  e)]  •  7  . 

If  solutions  y(x,  e)  of  an  equation  resembling  Eq,  (6.II) 


(6.11) 


[€L  +  €-“^x,  e)]  •  5 

are  known,  we  will  refer  to  Eq,  (6.12)  as  the  transformed  related 
equation,  for  it  is  readily  seen  that 


s 


y  =  s  •  y 


(6.15) 


satisfies  an  equation  of  the  form  of  Eq,  (6.4),  In  fact,  we  find 


R(x,  e) 

AAi- 


S  •  (J  -  U) 


(6.14) 
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which  has  the  desired  power  series  expansion,  since  the  exponential 
factors  in  the  S  ,  cancel  in  Eq.,  (6.l4)  while  the  (l  +  €  Q^)  matrices 

have  well-known  power  series  expansions  for  e  sufficiently  large.  Thus, 
we  regard  both  Eq,  (6,4)  and  Eq,  (6,12)  as  related  equations. 

Solution  of  the  Related  Equation 

Row  it  is  clear  that  the  matrix  ^  in  Eq.  (6.12)  is  to  he  chosen 
to  have  the  special  black-diagonal  form  of  the  ^(v  ^  l)  matrices; 
for  then  Eq,  (6,12)  decomposes  into  two  independent  differential 
equations 

^2  p 

+  [e  +  q(x,  e)  ]  y  =  0 

dx"^ 

y(x,  e)  =  y(x,  e)  •  e^ 
m-2 

q(x,  €)  =  y;  €-''i^3(x)  *  €-"^<x,  e)  (6.15) 


[eM(x,  €)  +  €"'^(x,  e)]  •  p 

m 


(6.16) 


where  "p  is  the  projection  of  y  onto  the  lower  (n  +  2) -space  and  where 
we  have  written 
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0 


0 


0 


J(x,  e) 


e(x,  e)  0 

0  Mfx,  e) 


Equation  (6.15)  is  the  type  discussed  by  Langer  and  McKelvey  and  we  will 
choose  e(x,  e)  so  that  it  is  precisely  their  related  equation.  We  can 
first  dispose  of  the  p  equation,  irrespective  of  the  nature  of  the 
turning  point,  i.e.,  irrespective  of  the  order  of  the  zero  of  Xq, 

Eq.  (6.7). 


Solution  of  the  (n  +  2) -order  system 


Now  Eq.  (6,l6)  is  of  the  type  considered  in  Sec.  5  and  4  except 
that  there  are  no  turning  points;  instead  the  eigenvalues  of  the  leading 
matrix,  n(x) ,  are  everywhere  identical,  and  series  solutions  can  be 
derived  by  the  method  described  by  Friedrichs,^  Thus,  we  set 


p(x,  e)  =  c(x,  e)  exp|e  J'  n(x')dx'J 


(6.17) 


whence 


dc 

dx 


=  Ni(x)  +  i  Np(x)  +  ...  +  — i- W  +iM]*c.  (6.18) 

*-  €  e  J 


We  now  demonstrate  that,  for  proper  choice  of  M,  Eq.  (6.l8)  has  a  complete 
set  of  n  +  2  independent  solutions  which  are  polynomials  in  l/e^ 

(6.19) 
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with  the  degree  p  as  yet  tinspecified.  Substitution  into  Eq_.  (6,l8)  yields 


(6.20) 


Equating  coefficients  of  €  ^  for  v  <  m,  we  obtain  relations  independent 
of  M, 


_ 

dx  ~  i-1 


SN  n*c  ,v<m 
>oM.+l  v-p’ 


(6.21) 


in  which  we  define  c,  =  0,  for  k  >  p,  the  degree  of  the  polynomial.  The 
first  of  these  relations  (for  v  ==  O)  is  seen  to  be  a  linear  homogeneous 
differential  equation  and  has  n  +  2  independent  solutions  which  form 
the  colvttnns  of  the  solution  matrix  ^(x),  which  can  be  specified  uniquely 
by  the  initial  data. 


(6.22) 


The  remainder  of  the  differential  equations  (6.2l),  i.e.,  v  ^  min 

(m  -  1,  p),  are  inhomogeneous,  with  inhomogeneous  teirnis  dependent  only 

on  c,  with  k  <  V.  Therefore,  these  can  be  solved  serially  in  terms 
k 

of  the  solution  matrix^  of  the  homogeneous  part.  In  particular,  we 


obtain  a  set  of  n  +  2  solutions  for  ea,ch  v,  namely 


(6.25) 


Now  it  can  be  seen  that  the  degree  p  of  the  poljmomial  c  cannot 
be  less  than  m  -  1,  for  Eq.  (6.2l)  would  yield  for  p<v  ^m-lan 
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algebraic  linear  relation  among  the  k  <  v,  which  is  not  in  general 
satisfied.  Hence ,  we  have  p  ^  m  -  1,  and  in  fact  it  suffices  to  choose 
p  =  m  -  1. 


How  if  we  remove  the  early  powers  in  l/e  from  Eq.  (6 .20)  through 
Eg.  (6.2l)  and  let  C  denote  the  matrix  of  polynomials  c,  viz,, 


(6.24) 


then  Eg.  (6.20)  is  found  to  determine  M 


M  = 


,-l 


m-J.  m-1 

g-v+1  ^  .  C  •  C“^ 

V=1  [I=V 


(6.25) 


Since  has  a  power  series  expansion  for  e  sufficiently  large ^  it 
follows  that  the  M  contribution  to  e)  has  the  power  series  expansion 

reguired  by  Eg.  (6,12),  A  complete  set  of  solutions  for  Eg,  (6.I8)  are, 
therefore,  known,  and  through  Eg.  (6.I7)  the  solutions  p  of  Eg,  (6,l6) 
"are  given. 


Solution  of  the  second-order  equation 


To  find  solutions  for  Eg,  (6.15)j  we  must  now  specialize  to  a 

particular  turning-point  problem.  Without  appreciable  increase  in  effort, 

2 

we  can  simultaneously  treat  both  types  of  extrema  in  c^ri.  In  order  to 
apply  the  formulas  of  McKelvey^*^  for  the  so-called  second-order 


80 


turning-point  problem^  we  transform  thus: 


d)  =  6€ 

6  =  i 

for  d^X^/dx^  >  0  at  X* 

6  =  1 

for  d^X^/dx^  <  0  at  X* 

Z  =  X  -  X* 


u(zj  co)  =  +  z,  ^/s) 

Po(“)  =  - 

Pj^(z)  =  -  6"\(x*  +  z) 

CO 

Q,(z,  a>)  =  -  q(x*  +  z,  &/(>)  =  ^  qj(z)co  ^ 


to  provide  a  notation  sljnilar  to  Ref.  10.  The  final  equation  defines 

the  low-order  terms  q.  in  terms  of  the  low-order  terns  of  Eq.  (6.15). 

D 

The  latter  equation  becomes 

.  [qSqCz)  +  &p^(z)  +  Q(z,  <£)]  u  =  0,  (6.27 

dz 

and  Eq.  (6.27)  satisfies  the  hypotheses  made  by  McKelvey.  In  particular, 
note  that  Pq(z)  is  non-negative  for  real  values  of  z.  [it  should  be 


^  The  reader  is  cautioned  that  several  specialized  symbols  are  introduced 
and  used,  particularly  between  this  point  and  the  end  of  this  subsection. 
Many  resemble  symbols  used  elsewhere  in  the  report,  but  differ  in  some 
detail.  Thus,  in  Eq.  (6.26),  u,  Pq,  and  p^^  are  not  to  be  confused  with 
mass  velocity  and  pressure.  In  addition,  the  independent jariable  z 
differs  from  that  of  Sec.  4  and  5  virtue  of  a  translation  of  the  origin 
by  but  the  symbol  is  identical. 
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recalled  that  we  restricted  attention  to  an  interval  containing  only  a 
single  zero  of  ^^(x).]  Following  McKelvey,  we  define 

<>(z,  0))  =  Pq(z)2 


(6.28) 


Note  that  <ti(z)  is  analytic  for  all  real  z  if  we  let  arg  4>  =  0  for 
z  >  0  and  arg  ^  =  n  for  z  <  0.  Thus,  arg  f  =  arg  a>  for  z  >  0  and 
arg  f  =  2n  +  arg  (d  for  z  <  0.  [For  arg(a)) ,  we  use,  of  course,  either  0 
or  Jt/2.]  Finally  we  note  that  II  is  regular  and  non-zero  at  z  =  0,  The 
introduction  of  the  variable  ^  has  the  character  of  a  "stretching  trans¬ 
formation"  since  it  appears  in  the  final  solution  below  as  the 
independent  variable  in  the  differential  equation  for  the  leading  part 
of  the  solution. 

Now  McKelvey  has  shown  that  the  behavior  of  the  solutions  of  Eq. 

(6,27)  for  large  |cd|  depends  essentially  only  on  Pq(z)  and  p^(z)  and  that 

the  solutions  are  asymptotically  given  by  the  confluent  hypergeometric 
17 

functions,  both  for  finite  z  and  for  the  neighborhood  of  z  =  0.  In 
particulax,  for  proper  choice  of  the  fiinction  e(x,  e) ,  Eq.  (6.15), 

McKelvey  has  shown  that  Eq.  (6,27)  has  series  solutions  of  the  form, 

u^(z,  (u)  =  T(z,  (c)"2[a(z,  a>)v^(z,  q)  +  B(z,  ^)v'(z,  a3)/a)]  (6.29a) 
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where  the  prime  denotes  the  derivative  with  respect  to  the  spatial 
coordinate  and  where  subscript  v  appears  explicitly,  not  only  to  label 
the  independent  solutions  but  also  to  show  that  the  polynomials  A  and  B, 


A(z,  od) 


m-2 


B(z,  (b) 


m-2 


5 


(6.29b) 


are  Independent  of  v.  The  leading  term  in  these  solutions  will  be  our 
principal  interest,  so  we  observe^^  that 


a^Cz)  =  cosh  j(z) 
bQ(z)  =  sinh  j(z) 


PiCzO  -  Pi(o)n^(o)/n\z 


(6.29c) 


It  is  important  to  note  that  j(z)  is  regular  at  the  turning  point, 
z  =  0,  We  have  for  T  the  determinant 


T(z,  cD)  = 


A  B/m 

a'  +  BS/ob  A  +  B'/tt) 


E  =  <ib%^(z)  +  ajp^(0)n^(0)/n^(z)  +  n'’‘^(z)/n(z)  +  k(a!)) 


(6.29d) 
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Fina3JLy,  the  functions  v^(z,  oi)  are  essentially  the  solutions  of  the 
leading  part  of  the  differential  equation  and  are  conveniently  written 

v^(z,  &)  =  k(cD)]  o  (6.29e) 


The  which  are  of  central  importance  in  the  sequel,  are  independent 

17 

solutions  of  the  confluent  hypergeometric  equation 


(6.50) 


for  the  special  value  p  =  . 

Now  the  choice  of  the  function  e(x,  e)  whereby  the  u^(z,  oj) 
become  the  exact  solutions  of  the  second-order  related  equation  need  not 
concern  us  (beyond  the  fact  of  its  existence) .  Instead  we  now  concentrate 
attention  on  describing  the  behavior  of  these  u^  functions  along  the  (real) 
z  axis,  particularly  with  regard  to  the  analytic  continuation  through  the 
turning  point  at  z  =  0, 

Now  the  neighborhood  of  z  =  0  is  seen  from  Eq,  (6.28)  to  correspond 

to  finite  values  of  ^  for  which  the  u^  functions  can  be  found  through  the 

l8 

power  series  expansions  of  the  confluent  hypergeometric  functions. 


k)  =  Mj^^p(t) 


M  (f)  =  1  +  f  + 

L  2p  +  1 


] 


P  = 


1 

4? 


1. 

4» 


(6.31) 
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Although  of  considerable  importance  in  proving  that  the  given  differential 
equation  and  the  related  equation  are  equivalent  asymptotically,  these 
relations  are  of  no  practical  importance  in  obtaining  the  desired  asymp¬ 
totic  estimates  for  the  u^  functions.  With  regard  to  the  latter  point, 
we  need  to  consider  only  non-zero  values  of  z  for  which  the  asymptotic 
behavior  of  the  W^  is  seen  to  be  important.  Thus,  for  large  |aj|,  we 
evaluate  Eq,  (6.29a),  using  for  T  the  result  from  Eq.  (6.29c)  and  (6.29d) 

T  ~  1, 

whence  we  obtain 

u^(z,  cc)  ~  n(z)  j^cosh  j(z)W^(t,kQ) 

dW  (t,k  )  n 

+  2  sinh  j(z)  - -  J  .  (6.52) 

Moreover,  the  derivative  of  u^,  which  enters  into  the  vector  solutions 
y,  Eq.  (6.12),  is  found,  on  using  the  differential  Eq.  (6.50)  to  evaluate 
the  second  derivative  of  W^,  to  satisfy  for  non-zero  z 

d-u  1  r  dW  (t,  k  ) 

^  ~  3)())nc^  I  2  cosh  o(z)  - - 

+  sinh  j(z)W^(t,  k^)  j  .  (6.55) 

Wow  rather  than  use  the  ^  hypergeometric  functions,  Eq.  (6.5l), 
it  is  more  convenient  for  purposes  of  asymptotic  analysis  to  use  the 

19 

Whittaker  functions  which  can  be  defined  as 
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(6.54) 


w. 


(f)  =  M.  (f)  +  — \  M 

^  rf4-^.k)  r(^P-k) 


r(l-p-k) 


r(|+p-k) 


and  particularly  the  va, riant  given  by  McKelvey 


10 


wjf)  =  W  ,  V  =  0,  ±  1,  ±  2,  ... 

V  Uil^  4 


0  =  (.1)''  (6-55) 

in  which  the  argument  of  the  product  of  complex  numbers  is  understood 
to  be  the  sum  of  their  arguments,  and  where 

arg(e'"^^^)  =  -  vjt# 

20 

These  have  the  asymptotic  form, 

+  0(l/t)],  for  f  on 

(v  -  3/2)1:  <  arg  f  <  (v  +  3/2)rt  .  (6.36) 

The  asymptotic  behavior  of  the  derivative  of  follows  from  the 
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derivative  of  the  Whittaker  functions 


dW, 


^  )  (k  -  |^)W  ^  -  [p  -  (k  -  I)  ]W. 


"df 


k,P 


Vl,p 


whence 


dW^ 

W 


for  t  on  Ey  . 


(6.37) 
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Now  it  is  evident  from  the  definition  of  the  intervals  that 
Eq.  (6.56)  and  (6.57)  are  not  sufficient  to  complete  the  expansions  of 
u^,  Eq,  (6.52)  and  (6.35),  for  z  both  positive  and  negative.  Therefore 

we  require  a  continuation  formula  for  the  functions.  For  this  puirpose, 

22  .  •  .  ■ 

we  employ  Kummer*s  "first"  formula  which  is  paxticularly  convenient  in 

the  form  given  by  McKelvey^^ 


\,p(«  = 


for  V  any  signed  integer  or  zero.  Using  this  result  in  conjiinction  with 

the  definitions,  Eq,  (6.34)  and  (6,55) >  and  employing  the  "reflection" 

25 

formula  for  the  gamma  functions,  we  obtain 


„  .  crkiti 

W  (f)  =  - ^ — - 

r(i-ak)r(3A-ak) 


,  2Gkrti„ 

+  e  W 


v+2 


W 


(6,38a) 


connecting  three  "adjacent"  functions  and 


o  •>  .  o  •  -bkrti 

=  -e"  +  /rA'  V'-;; — r  V5('^)  • 

r(i-ak)r(3/4-ak)  ^ 


(6.38b) 


It  is  evident  that  Eq,  (6.58a)  permits  the  continuation  of  the  asymptotic 
expression  for  from  to  the  intersection  n(S^_^^,  while  Eq, 

(6,38b)  permits  extension  to  (^(^^^2*  ^v+3^  *  sectors  are  found  to 

overlap,  and  further  continuation  is  given  by  application  of  the 
formula, 

(6.58c) 
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which  is  obtained  in  the  same  manner  indicated  above. 

Wow  it  can  be  shown  that  any  pair  of  four  "adjacent"  are 
linearly  Independent  [except  for  special  positive  values  of  k  =  -^  or  5/4 
(mod  l) ,  for  which  only  adjacent  pairs  are  independent].  It  is  convenient, 
for  the  purpose  of  finally  solving  the  stability  problem,  to  consider  the 
V  =  0  and  v  =  2  pair  and  to  obtain  from  these  u^  functions  the  correspond¬ 
ing  solutions  of  the  original  second-order  related  equation,  Eq.  (6.12), 

y  (x,  €)  = 


We  compute  now  asymptotic  expressions  for  y'‘'^'^on  z  >  0,  where  arg(f)  =  0 
or  n/2y  by  using  Eq.  (6.56)  and  (6.57)  directly  in  Eq.  (6.52)  and  (6.55) 
to  obtain 

z  >  0  (6.40a) 


where  ♦  has  been  evaluated  explicitly  from  Eq.  (6.7),  (6.26)  and  (6.28). 
For  z  <  0,  arg  ^  =  2jt  or  5^/2,  so  that  we  employ  the  continuation  formula 
(6,58b)  for  Wq  and  thence  expand  W^,  and  their  derivatives  to  obtain 


y'‘  '(x,  e) 


lU 


0 

0 

0 


Uy(z,  <n) 

€"^dUy/ dz 


0 

0 

0 


(6.59) 
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1 


y^°^(x,  e) 


-/cs/t 


lU 


0 

0 

0 


2ixie 


2kn:i 


r(i-k)r(5A-k) 


1 

Ks/tiu 

0 

0 

0 


In  similar  fashion,  we  obtain 


-*(2)/  \ 

(x,  e) 


1 

/Cs/t]U 

0 

0 

0 


-2rti 


r(i-k)r(5A-k) 


n(z)e2^ 


e) 


1 

.^Cs/TjU 

0 

0 

0 

1 

Ks/rju 

0 

0 

0 


z  <  0. 


(6.40b) 


(6.4la) 


(6.4lb) 


-*(l)  -*(3) 

¥e  could,  in  similar  manner,  determine  additional  solutions  y'  and  y' 
but  shall  not  do  so  here,  for  we  have  no  need  for  these  in  obtaining  the 
asymptotic  formula  for  6. 

We  have  now  obtained  a  total  of  n  +  4  solutions  of  the  related 
Eq.  (6.12),  and  we  denote  these  y^,  •••,  where 

e)  =  y^°\x,  e) 

y^ix,  e)  =  y^^\x,  e) ,  (6.42) 

while  the  remainder  are  formed  from  the  n  +  2  solutions  p(x,  e)  given  in 
Eq.  (6.17),  (6.19),  and  (6.25).  The  final  problem  is  to  prove  that  there 
exist  solutions  if  of  the  given  Eq.  (6.II)  which  agree  with  these  y^  up  to 
terms  of  order  e"™,  where  m  is  the  arbitrary  positive  integer  selected 
in  advance.  The  method  of  proof  follows  closely  that  of  Langer"^  and  is 
sketched  in  Appendix  D,  Based  on  this  proof,  we  have 

^y(x,  g)  =  yy(x,  e)  [1  +  0( €“"')]  (6.45) 

for  arbitrary  m. 

Asymptotic  Evaluation  of  6 

We  now  turn  to  our  main  task  of  obtaining  an  asymptotically  valid 

expression  for  ^  for  the  specialized  turning  point  under  consideration, 

2  .  .  . 

Again  we  return  to  the  specific  case  of  a  c^t]  profile  having  a  single 
maximum  (profile  M) .  It  is  seen,  then,  that  the  function  ^^(x),  Eq.  (6,7) j 
has  a  minimum,  whence  in  Eq.  (6,26)  and  the  following  discussion,  6  =  i. 
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Now  corresponding  to  each  of  our  solutions  Eq.  (6.45),  of  Eq.  (6.II) 
we  ha,ve,  hy  the  latter  equation,  a  solution  of  the  original  Eq.  (6.5). 
For  large  e,  it  is  seen  from  Eq,  (6,9),  (6.8),  and  (6,5)  that  only  the 
transformation  is  important ,  and  we  obtain 


(6.U4) 


for  all  X. 


To  find  9  asymptotically,  we  require  the  particular  Eq.  (6.44), 
which  remains  bounded  as  x  -+  We  contend  that  9^,  in  fact,  satisfies 
this  condition;  for,  writing  out  the  latter  explicitly,  using  ^  from 
Appendix  C  and  ~  y^°^  from  Eq.  (6.4o),  we  have 

^^-^[/%(z')dzj] 


^  _  z  k 


exp 


[ie(p^  +  [i-^)  +  ie1>(z0 

P3_(z0  -  Pi(o)n^(o)/n\z')^  ^  I 

2<D(zO  J  )  ’ 


(6.45) 


where  T  is  the  transform  matrix  Eq,.  (5*l)  and  where  the  variable  z  -  x  -  x* 
must  be  carefully  distinguished  from  x.  Since  for  x  >  x*,  ‘l’(z),  Eq.  (6.26) 
and  (6,7) ,  is  the  real  positive  quantity 
1  iM'. 


<l)(z)  =  p2  =  > 
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,  X  >  X* 


i6,k6) 


we  see  that  the  e  term  in  the  exponential  is 

-'X  /*X* 


/X 

p^dx'  +  ^ 


(6.47) 
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Thus,  6^  is  of  the  type  in  the  sense  of  Sec.  3o  Bringing  the  function 
of  X,  other  than  (T  •  e.),  into  the  exponential,  vre  obtain  the  e-independ- 
ent  part  of  the  latter  to  be 


/ 


;e  (x*+z')  -  ;e^(x*)n\o)'t>^(z')/  f  <P(z"')dz" 

_i - ± - ID - 1 - dz' 

2<t>(z') 


(6.48) 


where  the  point  z^  is  an  arbitrary  positive  number.  Now  by  Eq.  (6.29d) 
and  (6.26) 


kg  =  -  J  iij^(jr»)n’*(0) 


(6.I19) 


While  from  Appendix  C,  Eq.  (C.8),  (C.I9),  and  (C.20), 


(6.50) 


whence  it  follows  that  Eq,  (6,48)  has  as  its  x  derivative  as  in  the 

asymptotic  expression  (3,5)}  as  seen  in  Eq.  (3.6).  (This  is,  of  course, 
a  necessary  conclusion  since  the  present  asymptotic  expressions  must  agree 
with  those  of  Sec.  5  away  from  the  turning  point.)  Therefore,  the  present 
is,  outside  of  a  multiplicative  factor,  identical  with  0^  of  Sections  3 
and  4,  and  we  have,  as  in  Eq.  (4.6), 
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e(x,  e)  e^(x,  e). 


(6.51) 


but  now  without  restriction  on  x. 

The  solution  for  x  <  is,  therefore,  given  through  Eq.  (6.40b) 

to  be 


e) 


I- 


2rtie 


2fcn:i 


r(i-k)r(5A-k) 


(2ie)"^"^  f  (T  •  e  ) 

L  J 


+  e 


-Ukrti 


(6.52) 


X 

CieCii^+^g)  +  goldx' 

<t>(  z)  ^ 

Since,  for  x  <  x^,  ^(z)  is  the  negative  real  quantity, 

(^Cz)  =  Pq  =  ^^21  >  ^  ^  (6.55) 


the  first  contribution  to  Eq.  (6.52)  is  of  the  |J.^  type  and  the  second  of 

the  1J.2  type.  If  we  rescale  0  so  that  the  b^-type  contribution  is  simply 

T  •  e.,  for  X  =  0,  we  can  write  for  x  <  x* 

1 
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/X  ^ 

0 


/X 

(€^2  +  E22)<ix' 


p  >^-x*  n 

0!(e,  u)  =  -  /  *C2')dz  — ;- 


(6.51*) 

2k„  r(i-k„)r(3/4-k„) 

exp  I  (^>^X2  "'’  ^12^ 


exp  (6kQTti)  *0 


X^2  = 


x^(x*+z)  -  -6^(x*)n^(o)/n^(z) 
4>(z) 


This  is  precisely  of  the  form  of  Eq.  (4.24)  for  all  other  Class  II  values 
of  for  profile  M,  but  vrith  a  redetermination  of  a. 

Now  the  expression  for  V(tj  e)  (whose  zeros  determine  stability)  is 
to  be  obtained  from  Eq,  (6.54)  and  (6.45)  (the  latter  must  first  be 
rescaled  as  above,  of  course)  by  following  again  the  procedure  of  Sec.  5* 
It  follows  once  more  that  the  final  asymptotic  expression  is  given  by 
Eq.  (5.1?)  and  we  turn  now  to  consider  the  nature  of  Ot  at  the  current 
value  of  However,  there  is  no  need  to  locate  the  roots  of  L,  as  will 
be  seen  by  considering  the  magnitude  of  o:(e,  u) , 

Since  pure  imaginary,  |o:  j  ,  in  its  dependence  of  g,  is  seen 

from  Eq.  (6.54)  to  be  given  by  the  factor  One  finds  for  k^, 

Eq.  (6.^49),  using  Eq.  (6.28), 

k^  =  -  X^(x5^)/44''(0)  (6.55) 
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and  using  Eq.  (6.50)  and  the  values  of  and  Appendix  A,  this 

becomes 


^0  = 


^00  ”^^01 


_  2u^  ■ 

’  4  ♦'(0)tiu  (  L  ^ 


a*r 

qu 


(l-n)a.r^  dinq) 

Tim  ”  dx  » 

^x=x* 


(6.56) 


Since  is  positive,  Re(kQ)  is  clearly  negative,  provided  u  satisfies 

Re(u)  >  max  (0,  k^^/k^^) .  (6.57) 

For  such  t),  then,  0!  vanishes  in  the  limit  and  only  the  |J.^-type  solution 
contributes  to  6  asymptotically.  If,  however,  k^^  >  0,  then  in  the  range 

0  s  Re(o)  <  k^^/kQ^,  (6.58) 

0!  becomes  infinite.  In  the  expressions  for  6  and  L,  the  Pg-type  solution 

dominates  in  the  e  “  limit.  This  behavior  is  just  that  predicted  for 

profile  M  previously;  for  it  is  readily  seen  from  Eq.  (5 .l8)  that ,  if  k^ 

is  positive,  then  (in  the  Sec.  5  expression  for  a)  for  near 

i(c  Ti^)  is  large  and  positive.  Hence,  the  p„-type  contribution  also 
^  o '  'max  .  ^ 

dominates  for  just  below  the  special  value  of  Similarly,  if  k^  is 
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negative,  then  Pg  is  large  and  negative  and  the  p.^-type  contribution 
dominates.  In  view  of  the  criterion  for  instability  arrived  at  in 
Sec.  5,  viz.  Eq..  (5 .20),  the  criterion  for  instability  at  the  present 
special  value  of  is  the  inequality 


,^00 


(6.59) 


which  is  a  purely  algebraic  condition  on  the  equation  of  state  and  the 
chemical  kinetics,  as  opposed  to  the  integral  condition,  Eq.  (5 .20)  . 


Continuity  of  the  Asymptotic  Solution  for  Profile  M 

Having  arrived  now  at  a  stability  criterion  for  all  Class  III  values 

of  we  have  completed  the  treatment  of  profile  M,  except  that  we  still 

have  not  demonstrated  that  the  asymptotic  expression  for  V(t,  e)  behaves 

1. 

continuously  with  ^=€^  +  13  through  the  point  ^ 

of  the  continuity  of  V  requires  that  for  x  =  0  the  contribution 

should  vanish  when  ^  passes  from  Class  III  to  Class  Ij  for  we  recall 

from  Sec.  4  that  for  ^  in  Class  I  there  is  no  0^  contribution. 

From  the  discussion  of  the  magnitude  of  0£(e,  o) ,  it  is  eArident 

that  in  the  stable  case  (k^^  <  O) ,  the  contribution  does  vanish 

continuously  with  since  o:(e,  u) ,  Eq.  (5.18),  approaches  zero  at 

t  =  (c  Ti^')  .  In  the  unstable  case,  however,  it  is  clear  that  the 

=  i  'o'  'max 

contribution  vanishes  discontlnuously  with  ^  at  this  special  value 

i(c  r\^)  )  if  Re(u)  lies  in  the  range  Eq.  (6.58),  We  now  will 

o  rtiQiX 
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demonstrate  that,  for  less  "violent"  changes  of  t  = cC  +  h  (i.e.,  changes 

of  o  rather  than  0 ,  the  contribution  also  vanishes  continuously. 

Now  we  have  already  seen  that  by  increasing  Re(v))  through  the  value 

[see  Eq.  (6.58)],  the  desired  effect  takes  place.  Similarly 

for  Re(b)  <  k^^/kQ^,  we  will  now  show  that  the  magnitude  of  a  diminishes 

with  increasing  Im('o)  and  that  continuity  holds  in  this  direction  in  the 

T-plane.  To  study  the  magnitude  of  0!(e,  n) ,  Eq.  (6.54),  with  variations 

24 

of  Im('D),  we  first  employ  the  Legendre-Gauss  multiplication  theorem 
to  write 

r(i  -  Eq)  r(5/4  -  Eq)  =  r(i  -  . 

Consider  now  the  magnitude  of  Ci  for  both  e  and  Im(b)  =  large;  using 
the  above  identity  with  Eq.  (6.54),  we  obtain 

lo:(6,  n)  I  =  K^rd  -  2kQ)  exp  [2kQ  in  e  -  5kQ^jtP^]  (6.60) 

with  Kq  composed  of  various  factors  independent  of  and  e.  The  gamma 

25 

function  is  given  asymptotically  in  by  the  formula 

Ird  -  2  Re(kQ)  +  2ikQ^b^]  I  ~  exp  [-  2  Re(kQ)  in 

with  KLj^  representing  a  constant.  Combining  this  with  Eq,  (6.60),  we  see 
that  |0!  I  becomes  small  for  of  order  of  in  e,  i.e,, 

p.  >  P. 

1  1C 

p^^  =  in  e  .  (6.61) 
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The  continuity  of  the  asymptotic  theory  for  profile  M  is,  then, 
established, 

7.  A  -  B  DETONATIONS 


The  preceding  theory  is  here  applied  to  the  idealized,  one-reaction 

4 

A  -♦  B  detonations,  having  Arrhenius  rate  constant.  The  equations  for 
the  steady  flow  are  given  in  Ref,  4  and  are  repeated  here  only  in  so  far 
as  required.  The  notation  of  that  article  is  slightly  modified  here;  we 
note  particularly  that  in  the  present  notation  the  Mach  number  of  the 
shock  is  K  (rather  than  /c,  which  we  use  here  for  the  local  Mach  number) 
and  that  the  frozen  sound  speed  is  c^  (rather  than  c),  so  that  the  frozen 
sound  speed  in  the  quiescent  gas  is  c^_  (rather  than  , 

2 

Profiles  of  c  n 

o  ' 


One  readily  computes  this  quantity  in  terms  of  the  mole  fraction 
\  of  reactant  A  to  be 


2 

C  T1  = 
O  ' 


(1  - 


i^(x)  =  1  -  (1  -  k)  o/n 

03  =  {k^_  -  +  1) 

n  =  -  i)^/2(7o  -  1)^ 


(7.1) 
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where  Q  is  the  heat  of  reaction  in  units  of  RT_,  and  where  |  varies 
between  1  at  the  shock  to  some  equilihri'um  value  ^  0,  Since  X  decreases 
to  zero  monotonically  according  to  the  Arrhenius  rate  law 

dx/dt  =  -  fix  exp  (-E*/RT)  (7.2) 

(where  fi,  the  pre-exponential  factor,  and  E  ,  the  activation  energy,  are 

2 

constants) ,  the  derivative  of  c^t]  with  respect  to  x  has  the  sign  of 

(2oj|  -  1). 

A  particularly  simple  case  arises  if  the  ratio  of  (frozen) 
specific  heats  is  taken  greater  than  2,  for  then  (2cd|  -  l)  is  readily 
seen  to  be  negative  for  |  =  1  and  hence  for  all  0  ^  X  ^1,  irrespective 
of  detonation  velocity  and  irrespective  of  heat  of  reaction.  In  this 
case,  all  detonations  exhibit  profile  D  and  are  stable  asymptotically  in 
e. 

For  1  <  7  <  2,  we  examine  2cd|.  -  1  at  the  shock  (|  =  1)  and  find 

'  o 

it  to  be  positive  unless 

<  5/(2  -  7o). 

For  "normal"  values  of  y  ,  i.e.,  less  than  the  ideal-monatomic-gas  value 

'  o 

of  5/5,  it  is  not  expected  that  detonations  will  have  such  small  velocities, 
but  for  sufficiently  small  heats  of  reaction,  viz. 

Q  <  Q(7o) 

Q(7o)  =  7o(7o  +  1)/6(2  -  7o)(7o  ' 
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the  CJ  Mach  nttmber 


=  1  +  h  (1  +  +  ih^ 

h^  =  2(yl  -  1)q/7^  (7.4) 

is  less  than  5/(2  -  7^).  Thus  the  situation  for  Q  <  ^  is  that,  for 

p  p 

K  <  ^l{2  -  7^),  c^T)  decreases  with  distance  at  the  shock  and  throughout 

the  reaction  zone  and  exhibits  profile  D  once  more, 

2  2 

For  K  >  5/(2  -  7^),  c^Ti  is  increasing  at  the  shock,  but  not 

necessarily  throughout  the  reaction  zone.  Since  |  is  monotone-increasing 

in  X  (which  decreases  with  x)  and  since  (2d)|  -  l)  is  monotone-increasing 

2 

in  I,  it  is  seen  that  c^t]  has  at  most  a  single  maximum.  Now  the  value 
of  (2(o|  -  1)  at  equilibrium  is  negative  for  the  CJ  detonation  (for  | 
vanishes  at  \  =  0  in  this  case)  but  becomes  positive  for  sufficient 
overdrive  since,  by  Eq,  (7.I),  I  remains  unity  throughout  the  reaction 
zone  in  the  limit  K  .  There  exists,  therefore,  a  transitional 

Mach  number  K  ^  above  which  (2cjd|  -  l)  is  positive  throughout  the 
detonation.  One  finds 


2 

The  behavior  of  c  11  is  summarized  in  Table  III  for  all  values  of  the 
o  ‘ 

parameters.  It  should  be  remembered  that  the  possibility  of  instability 
is  limited  to  profiles  I  and  M. 


(7.5) 
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TABLE  III 


2 

Classification  of  the  profile  of  versus  distance  x  in  the  steady 
detonation^  according  the  types  in  Fig.  5>  A  -*  B  idealized 

detonation. 

Heat  capacity  Reduced  heat  of  Mach  number  of 

ratio,  7^  reaction,  Q  detonation,  /c_  Profile 


>  2 

<  2 


all 

<  ^(7q) 


all 


5/Ca  -  r^) 

5/(2  -  rJ 


2  2 
4  ^  < 


2  2  2 
'^-CJ  ^  ^  '^-t 


2  2 

K  ^  K 

-t  -t 


D 

D 

M 

I 

M 

I 


General  Considerations  of  the  Stability  Criterion 

The  implications  of  the  asymptotic  theory  can  be  found  by  the 
evaluation  of  exp  (^2)  ^1/^2  Class  III  values  of  for  specific 

values  of  the  reduced  heat  of  reaction  Q,  the  reduced  activation  energy 
Q*  =  E^/RT  ,  the  heat  capacity  ratio  7^,  and  the  degree  of  overdrive 
f  =  If  we  denote  the  value  of  |(x),  Eq.  (7.l)j  at  the  turning 

—  -Oo 

point  X*  by  |*,  then  the  correspondence  between  and  is 
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=  ■u^od(7^+  1)(1  -  (n|*)|*/(l  -  (7*6) 

Transformation  of  Eq.  (5.18)  to  the  variable 

i 

z(|)  =  [(1  -  a5|*)|*  -  (l  -  co|)|]^  (T»7) 

yields  for  Pg?  after  substitution  from  Eq.  (7.I)  and  (7.2), 

1  /-zCi) 

=  2co[(l  -  l2(z)dz 

I^Cz)  =  [(1  -  0)1)  (1  +  0)7^1)  r^(2a)|  -  1)“V2)  (7.8) 

/  t  \  2(1+0)7  1)  (1-0)1) 

K(2)  =  (r„  -  1)  -  1)  -  3J-  + 

|(z)  =  (20D)"^il  +  [1  -  ljO)(l  -  0)^*)l*  +  ikDZ^  ]2  I  . 


The  activation  energy  term  is  given  by 

+  (7  +l)^K^ 

p.g....g.q*ri-Kiy^g)(i-rt)r^  .  (7.9) 

(7o'^_+l) 

Now  it  is  evident  that  P2  becomes  large  and  positive  for  large  activation 
energy  and,  in  fact,  becomes  linear  in  Q,^,  For  any  for  which  is 
non-zero,  it  follows  that  the  detonation  becomes  unstable  for  all  suffi¬ 
ciently  large  activation  energies.  On  the  other  hand,  it  can  be  shown 
that  for  all  sufficiently  large  detonation  velocities 

|l^/L2  I  >  exp  (P2)  , 
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and  hence  detonations  are  stable  for  large  values  of  f. 

In  addition,  a  sufficient  condition  for  instability  can  be  derived 

1, 

by  a  comparison  of  the  derivatives  of  exp(P2)  and  L^/l*2 

The  derivatives  with  respect  to  itself  do  not  exist,  but  we  have 


1 — 

CVl 

2cd 

1  , 

O^ 

1 

mi). 

2(i)=o  " 

1 

e 

(7.10) 


where  z(l)  is  given  by  Eq,  (T*?)-  "the  former  exceeds  the  latter,  the 

exp  (^2)  curve  leaves  the  value  1  at  the  left  in  Fig*  7  above  the  L^/Lg 

curve,  whence  the  flow  is  unstable.  The  neutral  stability  point  with 

$  , 

respect  to  this  condition  yields  a  transitional  activation  energy 


for  each  Mach  niunber 


I"  (2a>.l)  (2x^-1) 


(7.11) 


which  is  independent  of  Q,.  For  >  Q^,  the  flow  is  unstable,  while,  for 

<  Q^>  exp  O2)  curve  begins  below  the  L^/L^  curve.  For  fixed 
values  of  Q.  a.nd  y  *  the  value  of  Q  below  which  the  latter  situation 
[with  regard  to  exp  (^2^  holds  for  all  detonation  velocities 

is  evidently  given  by  the  minimum  of  Q^.  Now  the  minimum  of  (for 
profiles  I  and  M)  occurs  at  the  smallest  detonation  velocity.  Thus,  for 
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Q  <  the  mnimuin  is  at  the  profile  D-M  transition^  where 

-  >'0^^  =  3(2’!r;7  »  ^  ^  ^  • 

and,  for  Q  >  Q,  at  the  CJ  detonation  velocity, 


(7.12a) 


In  order  that  a  system  with  given  Q  and  he  asymptotically  stable  for 
all  detonation  velocities,  it  is,  therefore,  necessary  that  Q  < 


Numerical  Results 


An  IBM-7094  computer  was  programmed  to  evaluate  Pg,  L^,  and  Lg 
throughout  the  turning-point  region 

(c  11^)  <  ^.  <  (c  ri^) 

'  o '  '+  1  '  o  •'max 


if.  2  2 

for  specified  values  of  the  parameters  Q,  Q  ,  f,  and  7^,  For  K.^  < 
the  value  of  k^,  Eq.  (6.56),  was  also  evaluated  for  L  =  (c  in 

order  to  ascertain  stability  at  the  maxim-um  of  profile  M.  A  plot  of 
exp  (Pg)  and  |L^/Lg  |  as  functions  of  as  discussed  in  Sec,  5  In 
connection  with  Fig.  7,  demonstrates  instability  by  the  presence  of  a 
region  for  which  the  exponential  lies  above  the  |L^/Lg  [  curve.  It  is 
to  be  noted,  however,  that,  in  the  absence  of  a  back  reaction  (B  -♦  A) , 
Pg(^i)  does  not  approach  -00  as  -*  but  attains  a  finite 

value,  as  seen  from  Eq.  (7.8).  This  value  of  (corresponding  to  the 
turning  point  being  at  x  =  “)  is  characterized  by  a  discontinuous 
behavior  of  0  and  V,  much  like  that  previously  discussed  for  profile  M, 


Inasimich  as  this  phenomenon  is  purely  a  result  of  the  reaction  models  there 

seems  to  be  no  compulsion  to  investigate  it  further* 

The  results  of  our  numerical  analysis  for  7^  =  1.2  can  best  be 

shown  in  the  form  of  the  neutral  stability  curves  in  the  (Q,  f) -plane. 

In  Fig.  8  are  shown  these  curves  for  the  activation  energies  of  Q  =50 

and  =  10.  It  is  seen  from  Eq..  (7«12)  that,  in  both  instances,  the 

2 

profile-D-to-profile-M  transition  locus,  viz.,  =  5.75 >  is  also  a 
neutral  stability  locus.  The  rightmost  neutral  stability  curve  has  a 
more  complex  structure  for  =  50.  The  lower  section,  i.e.,  the  portion 
below  Q  «  1,  corresponds  within  numerical  error  to  the  constant  Mach 
number  root  of  Eq.  (7.II)  for  =  50,  namely,  =  7*084.  In  the 
unstable  region  to  the  left  of  =  7.084,  the  exp  (pg)  curve  "begins" 
above  the  |l^/L2  1  curve,  but  for  some  below  the  latter, 

as  in  Fig.  7.  With  increasing  f,  found  to  decrease,  apparently 

merging  with  C-n  =  at  /c  =  7.084. 

II.  O  "T  « 

The  upper  section  of  the  rightmost  neutrally  stability  curve  does 

not  correspond  to  the  necessary  condition,  Eq.  (7.II).  For  fixed  Q 

above  Q  «  1,  the  exp  (Pg)  "begins"  below  the  |  Lj/li2  \  curve  w;hen 
P  1 

K  >  7.084,  but  for  the  range,  ^i  ^  above  the 

\\/\  I  curve.  These  are  then  unstable,  but,  with  increasing 

detonation  velocity,  the  unstable  range  of  decreases  [i.e., 

1. 

approaches  (c  t]^)  _  ],  finally  vanishing  (at  the  neutrally  stable  value 
of  f)  with 
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For  the  smaller  activation  energy,  Q,  =10,  the  unstable  regime 
is  considerably  smaller,  and  the  rightmost  neutral  stability  locus  lies 
to  the  right  of  the  root  of  Eq.  (7.II)  for  ^  =  10,  at  least  for  Q  s  0.5 
and  lies  very  near  to  this  root  for  Q  <  0.5.  If  a  "bulge"  of  the  type 
seen  for  Q,^  =  50  exists,  it  is  very  small  in  this  case.  It  should  also 
be  noted  in  connection  with  the  ot  =  10  results  that  the  previous  long- 
wavelength  calculations^  fo\md  instability  in  regions  of  asymptotic 
stability.  Thus,  the  Q  =  50  detonations  were  previously  found  to  be 
unstable,  apparently  for  all  degrees  of  overdrive  while  the  asymptotic 
theory  predicts  stability  for  all  f  >  1.04. 

The  most  astonishing  feature  of  Fig.  8  is  the  presence  of 

instability  for  a  range  of  f  values  at  all  positive  heats  of  reaction. 

4 

This  feature  was  not  present  in  the  long-wavelength  calculations  but 
its  correctness  is  indicated  by  the  fact  that  we  have  been  able  to 
calculate  an  unstable  root  for  Q  =  0.1,  f  =  4.0  by  means  of  the  numerical 
program  described  in  Ref.  4.  This  program  is  ineffective  for  large  values 
of  G  by  virtue  of  the  inability  of  the  program  to  integrate  the  0 
differential  equation  when  the  magnitude  of  the  eigenvalue  (having  negative 
real  part)  of  -  e)  5  Eq,  (2,4)^  is  large*  This  is  the  case  for 

large  €5  except  in  the  turning-point  region  of  the  T-plane.  Therefore , 
by  computing  and  as  well  as  P2?  L2  fo^  ^  particular  value 

of  for  which  the  condition  for  instability  holds ^  we  find  a  root  of 
the  asymptotic  expression  +  081j2  to  occur  at  the  pointy  [see  Eq*  (5  *18)] 
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Re(u)  =  - STTT — 

5'^i 

[  5rt/2P^(^.),  foTL^{±i^)>0 

e  =  ]  (7.15) 

(  n/2p.(q),  forL2(iCP<0 

Im(T)  =  eq  . 

For  i.  =  1.752818  u_^,  we  obtained  the  asymptotic  root 
e  =  15.77321 

T  =  0.001811  +  i  27.53626.  (7.l4) 

We  have  also  determined  an  "exact"  root  of  V(t,  e)  for  this  value  of  e 
by  detemnining  V,  by  means  of  the  numerical  calculation  of  Ref.  4,  for 
an  array  of  values  of  t  near  the  predicted  root.  Interpolation  yields 
the  root, 

T  =  0.0021583  +  i  27.556277.  (7.15) 

The  units  of  length  and  time  for  l/e  and  1/t,  respectively,  are  the  dis¬ 
tance  and  time  to  half-reaction,  as  in  Ref.  4.  The  accuracy  of  the  exact 
result  cannot  be  expected  to  exceed  the  third  digit  to  the  right  of  the 

decimal  point. 

The  transition  to  stability  for  all  Q  and  f  has  not  been 
intensively  investigatedj  but  for  <  4/5,  it  follows  from  Eq.  (7.12) 
that  the  sufficient  condition  for  instability  is  not  satisfied, 
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irrespective  of  Q  and  f.  Numerical  evidence  suggests  that,  as  the 
activation  energy  approaches  Q,^  =  4/3  from  above,  the  rightmost  neutral 
stability  curve  approaches  the  profile  D-M  transition  locus.  Hence 
for  <  4/5  it  is  believed  that  all  detonations  in  this  system. are 
stable  for  small  transverse  wavelengths.  No  attempt  to  prove  this 
conjecture  has  been  made,  but  it  has  been  proved  to  be  at  least 
approximately  true. 

One  final  point  seems  worth  mentioning,  namely  the  approach  to 
0,  =  0,  the  step  shock.  The  exact  expression  for  V(t,  e)  is  known  for 
the  step  shock,  of  course,  namely 

V=eL^(C), 

and  this  has  no  roots  on  the  imaginary  ^  axis.  As  described  above  in 
conjunction  with  Fig.  8,  the  asymptotic  theory  fomally  predicts  for 
Q  =  0  a  range  of  unstable  shock  velocities,  dependent  upon  the 
■activation  energy,  >  4/5,  In  the  Q  =  0  limit,  however,  it  can 
readily  be  shown  that  -*  0  and  ->  Thus,  according  to  Eq.  (7. 13), 
the  rate  of  growth  of  the  predicted  "instability"  is  zero  while  the 
wavelength  vanishes.  Thus,  while  the  asymptotic  theory  is  not  exactly 
continuous  at  the  Q  =  0  point,  it  does  not  predict  an  incorrect  result, 

8.  DISCUSSION 

It  is  clear  from  the  numerical  results  for  the  idealized  system 
that  the  asymptotic  theory  provides  considerable  additional  information 
on  the  behavior  of  detonations  and.  obviously  should  be  an  integral  part 
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of  any  stability  calculation.  In  addition,  the  possibility  that  rather 
small-scale  transverse  disturbances  [the  wavelength  of  the  root,  Eq, 
(7.15),  is  about  0.1  in  units  of  the  distance  to  95fo  complete  reaction] 
are  unstable,  even  though  long  wavelengths  are  stable,  is  certainly  not 
contrary  to  experimental  evidence.  The  latter  seems  to  indicate  the 
presence  of  small-scale,  transverse  inhomogeneities  even  when  the 
conditions  are  well  removed  from  those  associated  with  spin. 

Although  the  present  theory  does  not  completely  describe  the 
stability  behavior  of  a  detonation,  it  does,  nonetheless,  provide  a 
computationally  simple  criterion ,  the  fulfillment  of  which  is  a 
sufficient  condition  for  the  growth  of  disturbances.  The  application  of 
the  present  theory  would  appear  to  be  entirely  feasible,  irrespective  of 
the  number  of  chemical  reactions  and  the  complexity  of  the  equation  of 
state;  in  this  respect,  it  might  well  represent  the  only  theoretical 
test  of  stability  which  is  available. 

In  connection  with  two-  or  three-dimensional  calculations  of  time-_^ 
dependent,  reactive  flows  (should  such  ever  become  feasible) ,  the 
presence  of  asymptotic  instability  has  the  greatest  importance.  In  a 
calculation  involving  no  transport  effects,  continued  refinement  of  the 
mesh  would  evidently  increase  the  number  of  permitted  unstable  modes, 
and  a  mesh  extrapolation  would  apparently  be  meaningless.  A  realistic 
inclusion  of  viscosity  presumably  would  be  required. 
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APEEMDIX  A 


The  matrices  A  ,  A  ,  and  B  are  given  in  Ref.  5  in  terms  of  functions 
of  the  state  variables  in  the  steady  flow.  The  Inverse  of  A  is  readily 
found  to  be 


(A.l) 


where  I  denotes  the  unit  matrix  of  order  n,  m  is  the  mass  flux  u/v,  and 
/<-n 

p  is  the  pressure.  Subscripts  S,  v,  and  \  denote  partial  derivatives,  with 
the  remaining  variables  of  the  set  held  fixed. 

The  calculation  of  Eq,.  (2. 76),  involves  a  matrix  product 
yielding 
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(A.2) 


The  eigenvalue-eigenvector  problem, 

=  0 


(A.J) 


proceeds  by  finding  the  determinant  of  the  coefficient  matrix  to  be 


1#! 

T  u 


(A.4) 


This  has  roots  as  follows; 


s] 


^5  "  ^4  “  =  V4  " 


(A.5) 


s  =  [^^  +  cl^r  . 
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The  determination  of  the  corresponding  eigenvectors  is  somewhat 
laborious  but  straightforward.  The  T^  form  the  columns  of  the  matrix 
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The  calciaation  of  the  matrix^,  Eq.  (5.2)  is  a  laborious  procedure, 
■which  we  now  indicate*  We  had 


>  =  [A"^  •  (B  +  hi)]' 

»1  fmX  ^ 


(A.8) 


whence  we  write 

E  =  E^°^  +  hE^^^ 


E 


(o)  ^  .  [w 


e(1)  .  (a“^)  '  •  T. 

/w^^  jMk. 


(A.9) 


The  matrix  W,  is  obtained  from  Eq.  (A.l),  and  matrix  B  is  given  in  Ref.  3 
in  terms  of  the  gradients  in  the  steady  flow  and  certain  rate  coefficients. 
The  gradients  can  be  written  as  products  of  thermodynamic  functions  and 
the  chemical  rate,  so  that  B  can  conveniently  be  separated  into  a  term 
linear  in  the  rates  and  a  term  linear  in  the  thermodynamic  derivatives  of 
the  rates.  Thus  we  write 


(A.IO) 


with  B_  vanishing  at  equilibrium.  Explicitly 


and  has  hut  two  non-vanishing  rows, 
ywvL 


\  1  =  1.  2.  5 


•  B,  =  [r  ,  0,  0,-  r„,  rr*]  .  yp- 


(A.12) 


-  Ir,  Oj  0» 


Here  a  is  the  important  thermodynamic  function  defined  in  Ref.  7 j  ^  is 
the  gradient  operator  in  X-space^  also  given  in  Ref»  7?  "while  subscript  X 
denotes  the  partial  themodynainlc  derivatives 5  as  previously  noted*  In 
addition^  we  have  the  following  identifications:  T  is  the  thermodynamic 
temperature,  7^  is  the  usual  ratio  of  frozen  specific  heats,  AF  is  the 
free  energy  increment,  and  is  the  frozen  expansion  coefficient 

Po  =  (^/5T)p,>: 
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Miltiplication  by  A”^,  Eq.  (A.l)  with  subsequent  transposition  yields  a 
corresponding  pair  of  terms, 


ill 

=  w  +  w 

/vLlO  Zll 

/JilO 

1 — 1 

I — 1 

II 

(A.15) 

=  [a"^  •  B- ] '  . 

A^ll 

A/v5{  AmI 

explicit  formulas  for  each  are 

w  =  i 

AviO  U 


r  2(1. ri)  ^sl  a*T  ^  r»VT^  ^2-1)  ^sl  ma-r  ^  mr*^ 

■  — J  — "  - ^ J  T“  ^  ^ 


(2~3il)^*r 
mT] 

0 

.■;  (l-n)r-9(p  ) 


VOT 

“  ^ 

'  VO  Tim 


2(1-ti)  OT 

- 2 - 


ucj.r  ' 


Tjin 


^  m  vaT 
°  ° 


^  T 

uT 


7  vcr*r 


(*)x-']  -s"['-%>--s(r)j-]  »  (r)j-  » 


(A.14) 


(where  is  the  frozen,  constant  volume  heat  capacity);  and^^^  has 
three  non-vanishing  rows. 


Il6 


i  =  2,  5 


W 

Zil 


A 


w 

<illl 


(A.15) 


The  calculation  of  Eq.  (A.9) ,  involves  differentiation  of  £^as 
well  as  matrix  multiplications,  and  we  display  only  the  final  result  for 
several  elements.  Thus,  we  obtain 

^  [_  T  V  J 


11 


1  .  2-ti  ^ 
4ti(1-ti)  dx 


(.  r  j  TP^o*r  a«r  "I 

ilL.  _i_  ^  -  —12 _ I  - 

2tis  LI-ti  dx  ”  Tu  "  ni  J 


1  d  jjn  s 

2  dx 


/c^+s  Ccr*r 
*  i^+Ks  squ 


(A.16) 


,(1)  =  .  < 


11 


uns 


For  Egg  one  obtains  a  result  equivalent  to  "when  the  square  root  s 

is  replaced  by  -s.  The  off— diagonal  elements  of  not  enter  into  our 

final  discussion  in  any  important  way  but  it  is  of  interest  to  note  that 
the  element 
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E, 


(0)  _ 


21 


a*r 


2tiu 


(2  -  Ti) 


(l-Tl)TlVPg]  2-ti  dn 

T  "  iJ-Ti(l-ri  )  dx 


is  singular  at  a  Class  II  turning  point,  where  -  K-s")  vanishes.  It  is 
seen  from  Eq.  (3.7)  that  the  first  correction  term  to  is  singular, 
then,  even  though  the  leading  term  is  regular. 
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APPENDIX  B 


In  this  appendix,  we  prove  that  unstable  roots  of  V(t,  e)  occur 
in  conjugate  pairs  in  the  T-plane.  For  this  purpose,  it  is  sufficient 
to  show  that,  by  suitable  choice  of  the  arbitrary  multiplicative  constant 
in  0(x,  T,  g),  V(t,  g)  is  an  analytic  function  of  t  and  is  real  for  all 
real  values  of  t. 

Now  V  is  given  by  Eq,  (2 .5)  which  can  be  written  in  the  following 
fashion,  using  the  vectors  g^,  g^,  h^,  and  h^  of  Ref.  5: 


(B.l) 


From  Appendix  A,  Eq,  (A.l) ,  it  is  seen  that  the  products  of  A“^  with  the 
g  vectors  have  the  same  zero  elements  as  the  g's  themselves.  It  follows, 
then,  that  V(t,  g)  is  real  and  analytic  (for  real  t)  if  0  is  analytic  and 
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Re[^(x,  T,  e)]  •  =  0 

Im[^(x,  T,  e)]  •  =  0,  i  =  1,  2,  4,  5, 


(B.2) 


Now  the  differential  equation  for  0,  Eq,  (2,4),  has  the  coefficient 
matrix  -P^(x,  t,  e) 


P'(x,  T,  €)  =  .  [tI.  UA,.  B])' 


(B.3) 


which  has  real  part 


{a"^  •  [tI  +  B]}' 

fNt 


(B.4) 


and  imaginary  part 

€(A“^  •  A  ) ' 

^y 


(B.5) 


From  the  fact  that  B  has  only  zero  elements  in  column  5  and  in  row  5 
[See  Appendix  A,  Eq,  (A,ll),  (A.12)],  it  follows  that  the  real  part  of 
is  diagonal- with  respect  to  element  5>  i.e,,  it  has  the  same  block 
diagonal  form  of  , 


:Mi2 

0 

\4 

“15 

■^2 

0 

^4 

0 

0 

“33 

0 

0 

\l 

0 

M44 

“A5 

^2 

0 

^4 

•^5 

(B.6) 


On  the  other  hand,  the  imaginary  part  is  foiond  to  be  of  just  the 
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complementary  form,  i.e.,  it  has  non-zero  elements  only  in  the  off- 

diagonal  positions  of  row  5  and  column  5.  Hence,  if  0  has  the  form  of 

Eq.  (B.2) ,  then  d^/dx  also  does,  and,  hy  repeated  differentiation,  so  do 

a.l  1  its  derivatives.  If,  therefore,  0  has  the  form  of  Eq,  (B,2)  for  any 

X  ,  it  does  so  for  all  x, 
o’ 

Unfortiinately  the  only  point  in  the  reaction  zone  where  0  can  he 
determined  is  in  the  equilibrium  state,  x  -*  <».  It  turns  out,  however, 
that  the  property  (B.2)  holds  in  this  limit  and  its  validity  there  is 

sufficient  to  establish  it  everywhere.  We  begin  by  demonstrating  that 

3 

the  asymptotic  form  of  0,  as  x 

?(x,  T,  e)  ~  v(t,  €)  (B.7) 

has  the  property  (B.2).  Now  IJ.(t,g)  is  the  eigenvalue  of  -^'(“,t,  e) 
having  negative  real  p9«rt  and  v(t^  e)  is  the  associated  eigenvector* 

Since  P  is  analytic  in  ji  is  also  analytic  and  v  can  be  chosen  to  be 
analytic. Now  the  complex  conjugate  of  has  a  single  eigenvalue  with 
negative  real  part*  It  can  be  seen  from  the  properties  previously 
ascribed  to  its  real  and  imaginary  parts^  Eq.  (B.4)  and  (B.  5),  that  the 
following  relationship  holds: 

-  T,  e)  •  v^(t,  e) 

''l 
^2 

v^(t,  g)  =  - 
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■where  *  denotes  the  complex  conjugate.  Thus 

^*(t,  e)  =  h(t,  €)  (B.9) 

and  |i  is  real.  Moreover,  v*  is  also  an  eigenvector  of  -P%  whence  the 

1  AN 

eigenvector 

VgCT,  e)  =  i[v  +  V*  ■  -  ^f)l  (B.IO) 

has  the  property  (B.2), 


To  complete  the  proof,  we  obtain  an  explicit  solution  for  0, 
valid  in  the  neighborhood  of  equilibriim.  To  do  this,  define  a  new 
independent  variable 


(B.ll) 


where  S(x)  is  the  entropy  in  the  steady  reaction  zone  which  is  constrained 
to  increase  monotonically  with  x  by  virtue  of  the  second  law  of  thermo- 
dynamics.  Thus  z(x)  varies  monotonically  between  1  at  the  shock  and  0 
at  X  =  00. 


The  differential  equation  for  ^  becomes 


T,  e) 


n  =  P'(x,  T,  g)  (B.12) 

where  $  is  the  entropy  production,  -  £^*t/Tj  a  non-negative  thermodynamic 
function.  The  "z  =  0  point  is  a  singular  point  of  Q  by  virtue  of  the 
fact  the  $  vanishes. 
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fr 


Now  throughout  the  reaction  zone^  the  therrnodynarnic  functions  of 
state 5  under  the  steady-state  constraints,  were  shown  by  Wood  and 
Salsburg"^  to  be  regular  functions  of  the  progress  variables  so  that 
P  itself  is  regular  at  z  =  0.  Near  equilibriinn,  the  entropy  production 
$  is  of  second  order  in  \  -  \(°°) ,  and  the  entropy  gradient  in  \-space  is 
of  first  order.  Thus  \  -  X(“>)  vanishes  linearly  in  and  Q  has  a 
simple  pole  at  "z  =  0.  In  some  neighborhood  of  “z  =  0,  the  series 


n  = 


(2)'^ 

k=0  ^ 


(B.13) 


converges,  and  the  leading  term^  is  proportional  to^'(“,  T,  e). 

Convergent  power  series  solutions  of  Eq,  (B,12)  can  be  derived  by 

27 

application  of  the  theory  of  singularities  of  the  first  kind  ,  and  the 
lone  solution  bounded  at  z  =  0  can  be  picked  out,  viz. 


■t(z,  T,  e)  =  z^ 

Here  we  have  employed  fi  and  Vg  for  the  eigenvalue  and  eigenvector  of 
so  that  f  is  real  and  positive  [opposite  in  sign  to  |J.  of  Eq.  (B,7)]  and 
Vg  is  given  by  Eq.  (B.IO)  to  have  property  (B.2) .  Now  by  Eq.  (B.l4)  the 
quantity 

@  =  (z)"*^  ?(z,  T,  €)  (B.15) 


I  +  z 


k^O 


r.  -k 

Si? 


is  regular  at  z  =  0,  and  its  derivative 


( 


n 


(B.l6) 
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will  also  have  the  property  (B.2)  at  this  point,  as  follows  from 
the  previously  discussed  properties  of  P^,  In  fact,  all  the  derivatives 
obviously  do^  and,  since  these  are  precisely  the  higher  coefficients 
in  Eq.,  (B,14),  it  follows  that  0  itself  has  the  desired  property. 
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APPEKDIX  C 


The  transformations  of  Sec.  6  begin  with  ■which  satisfies  Eq^.  (6,6) 
[jo  -  *  (I'l  + 

The  simplest  route  to  finding^  is  by  first  diagonalizing  the  coefficient 
matrix  on  the  left.  We  have,  from  Appendix  A, 

T“^  •  -  -1(1^,  +  ^p)!]  *  ^i^)»  •••> 

(C.2) 


and  we  require  the  transform 

Diag  [-  ^A2l/2>  •••>  n(x)]  • 


il0‘ 


(C.5) 


One  readily  finds 


1 

2 


tiu/2Ks 

iTu/2k:s 


Thus 

0^  =  T  •  H 

iSO  I"*  /«*' 

has  columns 


^  * 


in+2 


0 

mr|/(l-Ti) 

- 

0 

11 

1 — I 
<(U 

i 

’  A  •  “ 

0 

0 

-  nPg/m  - 

0 

-  _ 

e .  = 

3  ‘ 

T  -  e.  , 
3  ’ 

j  =  5,  ***,  n 

+  4. 

(c.^) 


(C.5) 


(C.6) 
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Moreover j  since  we  had 

E  =  T"^  •  [0^  .  T  -  dT/dx] 

MVi  AW  MLt  AW 

we  have  for  the  matrix  H,  Eq.  (6,6) 
H 


(C.7) 


Thus  H  differs  from  E  only  in  the  first  two  columns. 


H 


A 


2  +  ^12  +  ^21  ^22^ 

-  ~  ^21  ■  ^22^ 
“  2'  ^  ^32^ 

"  ^  (^41  ^42) 

-  i  (^51  ^  ^52) 


H 


A 


qu 

t'^U. 

"^12'^  ^21  “  ^22^ 

i  (E 

h  ■ 

Ei2  -  E^^  +  E22)  -  ^  An(ilu/Ks)/dx 

qu 

("31 

“  %2^ 

qu 

“  V 

qu 

27^ 

'®51 

-  E^g). 

(0.8) 


It  should  be  emphasized  that,  despite  the  appearance  of  JE  and  s  in  Eq. 
(C.T)  and  (C.8),  H  is  a  regular  function  of  x  by  virtue  of  Eq.  (6.6), 
(C.6),  and  (A,6) , 

In  order  to  prove  the  existence  of  the  transforms  Eq.  (6.IO) ,  we 
employ  induction.  Thus  we  first  show  that  the  transformation 
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♦(0)  =  S  •  ♦(!) 


Si  =  [I  +  exp  g^dx' 


yields 


4F-  =  ^  ^  e’V^^Cx,  €)]  .  ? 


Wow  differentiation  of  Eq,  (C.9)  yields 


(1) 


(C.9) 


(c.io) 


dilf 


-  ^4.  *  ^(1) 


(C.ll) 


-1^  N  -lx 


Using  the  identity 

-1^  N-1 


(J  *  '  =  i -  '  Si  *  Sl  •  (i "  ’ 


we  have 


-  ^0  ^  ^  *  ^0  S' " 

•  ^  -  ^  *  (i  ^  So  • 

-  +  (^  +  e  ^)  •  (^  •  ^-d^/dx)]  .  (C.12) 

The  terms  within  the  bracket  have  the  desired  power  series  expansion 
irrespective  of  the  choice  of  and  so  that  we  need  only  show  that 


So  *Si"/2q  *So’’’S"  %£  Si 


(C.15) 
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for  suitable  choices  f or  ^  and  g^. 

To  solve  Eq,  (C,15),  it  is  convenient  to  drop  the  subscripts  on 
and  write  the  matrices  in  block  form, 


wll 

Sl2' 

1 - 

0 

0  H 

H  = 

&2_ 

11 

1^2. 

where  is  a  2-by-2  matrix  +  2).  Then 


Eq.  (C.15)  becomes 

&1  -  -  &1  •  £11 "  Jii  •  Sii  =  iii 

-ai2  -162  =  0 

&l -&1 

i?22  ■  ®oJn+2  "-Sea  '-£22  *£22  '  £22  '£22  •  (c.l4(l) 

Since  ^22  is  n(x)^^^2>  'w®  have 

Ssi2  *  -^22  ^  ^^^^^>812 » 

whence  Eq,  (C,l4b)  yields 

=  -  [ifn  -  •  &2  •  «=-i5) 


The  inverse  matrix  on  the  right  exists  by  virtue  of  the  fact  that  the 
eigenvalues  of  L°  are  distinct  from  n(x) ,  Similarly,  Eq.  (C.l4c)  has 
solution 


128 


(C.l6) 


•  fin  -  • 

In  view  of  the  character  of  evident  that  ^2 

disappears  from  Eq.  (C.lM)  and,  thus,  does  not  appear  in  Eq.  (C.l4)  at  all. 
Hence,  we  can  simply  write 

*822  ""  ° 

J^2  “  ^0mi+2*  (C»1T) 

Finally,  Eq.  (C.l4a)  becomes  in  component  notation  (i.e.,  dropping 


subscripts  from  and^^,  but  inserting  subscripts  for  the 

individual  elements) , 

four 

-  go  +  ^0^2  +  ^1  =  ° 

(C.l8a) 

0 

11 

OJ 

+ 

(C.l8b) 

^1  '^0*^2  “  •^o'^l  “  "  \ 

(C.l8c) 

H22  -  §0  "  ^1  "  ^0^2  " 

(C.l8d) 

Combining  Eq.  (C.l8b)  and  (C.l8c),  we  readily  obtain 

"  ^0^  "  ^1 

S.1"  ^2  "  ^  * 

(C.I9) 

Adding  Eq.  (C.l8a)  and  (C.l8d),  we  get 

go  =  i  (Hii  +  H22) 

(C.20) 

which  serves  to  determine  g^.  The  proof  is  completed  by  choosing 
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(C.21) 


-  1 

^  "  ‘^0  ^ 

%.1  =  ^2 

Q22  =  0  • 

To  complete  the  inductive  proof  of  Eq,  (6.IO),  we  assimie  its 
validity  for  k  -  1,  with  k  S:  2,  and  examine  the  next  transformation 

^(k-l)  ^/^k  ’  ^(k) 

dx  ^  ’  ^(k)  (C.22) 

lik  -  (£*  ■  t^4-i  •  'i" 

Employing  the  identity 

=  i-  •  <v  ’ 

we  have,  from  the  properties  of  ^  Eq.  (6,10), 

4  =  -  &  •  io  -  io  •  i  -  ViS 


We  do  not  write  out  the  terms  in  the  bracket,  but  it  is  evident  from 

Eq.  (C.22)  that  this  remainder  has  a  power  series  expansion.  The 
-k+1 

coefficient  of  e  ,  on  the  other  hand,  agrees  with  that  for  k  =  1, 
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Eq.,  (0,12)5  except  in  the  subscripts  on  the  unknowns  and  the  symbol  for 

fk  1) 

the  "knowis"  H  and  UX  whence 

-  ik  •  io  *  io  •  Jk  -  ik 

has  a  solution  of  precisely  the  same  form  as  before.  The  proof  is 
therefore  complete. 
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APEEiroiK  D 


The  proof  of  the  asymptotic  equivalence  of  the  (transformed) 
given  differential  Eq.  (6.11)  and  the  (transformed)  related  Eq.  (6.12), 
l.e.,the  proof  of  Eq.  (6.43),  Is  briefly  outlined  here.  It  depends 
upon  consideration  of  the  Integral  equation  for  the  solutions  of  Eq. 

(6.11), 

(x,  e)  =  y^(x,  e)  +  6““y(x,  e)  •  f  Y"^(s,  e)  •  r(s,  e)  •  ^^(s,  €)ds 
r(x,  e)  =  U(x,  e)  -  J(x,  e)  (D.l) 

aa- 

where  Y  denotes  the  matrix  of  coliimns  y. ,  Eq.  (6.42),  and  where  the 
lower  limit  of  the  Integral  on  the  right  may  depend  upon  the  "matrix 
element"  within  the  integral;  i.e.,  if  the  integrand  Is  written  as  a 
sum,  the  lower  limit  may  vary  from  term  to  term.  The  validity  of  Eq. 

(D.l)  is  readily  established  by  differentiation. 

.  In  order  to  establish  Eq.  (6.43),  it  is  conducive  to  clarity  to 
consider  the  specific  case  6  =  i  in  Eq.  (6.26),  corresponding  to  a  turning 
point  X*  at  which  c^q  is  a  maximum.  In  fact,  we  will  simply  restrict 
our  attention  to  profile  M,  Fig.  5j  which  there  are  no  other  turning 
points  on  the  x  axis.  By  virtue  of  this  specification,  the  "stretched" 
variable  t,  Eq.  (6.28)  has  argument  it/2  for  z  >  0  and  5it/2  for  z  <  0. 
Turning  to  the  asymptotic  expressions  Eq.  (6.4o)  and  (6.4l),  we  find  that 
the  behavior  of  these  functions  with  increasing  e  (i.e.,  increasing  |a)|) 
depends  on  the  terms  rather  than  on  the  exponentials  in 
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To  prove  the  desired  asymptotic  result,  it  is  important  that  the 
set  of  n  +  4  known  solutions  of  the  related  equation, 
y^(i  =  Ij  ••.,  n  +  4),  be  "asymptotically  distinguishable"  for  all  x. 

This  property  is  illustrated  by  the  solutions  y^  and  yg  if  we  assume 
that  Re(k)  <  0,  From  Eq.  (6.4o)  and  (6,4l),  it  follows  that  y^  =  y^*^^ 
is  recessive  relative  to  for  z  >  0,  while  for  z  <  0,  y^^^ 

is  recessive.  By  virtue  of  the  fact  that  one  is  recessive  while  one  is 
dominant  for  all  z,  these  are  said  to  be  asymptotically  distinguishable. 
For  the  case  of  Re(k)  >  0,  y^^^  and  y^^^  remain  linearly  independ¬ 
ent  solutions.  However,  the  leading  part  (for  large  e)  of  each  is  of  the 
dominant  type,  as  is  seen  from  Eq,  (6,4o)  and  (6,4l),  Therefore, 
asymptotically  they  are  not  distinguishable  and  are  not  suitable  in  the 
proof  of  Eq,  (6,45),  In  this  instance,  it  can  be  demonstrated  that  the 
pair  y^^^and  y^^^  given  through  Eq,  (6,59)  distinguishable  and  should 
be  used  for  y^  and  y^  in  the  proof,  1,6,,  in  place  of  Eq,  (6,42), 

For  y^  and  y^  in  the  case  Re(k)  =  0  and  for  the  solutions  of  the 
(n  +  2) -order  related  equation,  y^,  ••*,  "the  concept  of  asymptotic 

distingul shability  is  not  particularly  important J  for  the  various 
solutions  do  not  become  recessive  or  dominant  for  large  e.  Instead,  the 
leading  parts  of  these  solutions  in  the  €  -»  “  limit  remain  linearly 
independent , 

In  view  of  the  distingulshability  of  y^  and  as  defined  in 
Eq,  (6,42)  for  Re(k)  <  0,  we  will  restrict  our  attention  to  this  case  in 
the  remainder  of  this  Appendix,  The  Re(k)  >  0  case  is  equivalent,  once 
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and  have  been  redefined  in  the  manner  indicated  above  o 

The  proof  is  based  upon  an  adaptation  of  the  theorem  used  by 
9 

Langerr  Given  the  integral  equation 

/X 

K(x,  s,  e)  •  f(s,  6)ds  (D.2a) 

with  the  order  relations 


g(x,  e)  =  h(x,  6)0(1) 


then 


(D.2b) 


f(x,  e)  =  g(x,  e)  +  h(x,  e)0(€“"^). 


(D.2c) 


The  norm  to  be  used  for  the  arrays  is  the  maximum  absolute  value  of  the 
elements.  The  proof  of  the  theorem  is  essentially  given  by  hanger. 

To  begin  with,  then,  we  require  the  kernel  of  our  integral 
equation  in  somewhat  greater  detail.  In  particular,  the  inverse  of  Y  is 
required.  From  Eq,  (6,17),  (6,24),  (6,39) >  a^d  (6,42)  we  obtain 


0 


u, 


0 


2^^(x,  e)  exp  -  en(x')dx' 


(D.5) 


where  F  is  the  determinant  of  the  upper  2-by-2  part  of  Y  and  hence  is  e”' 
times  the  Wronskian  of  the  solutions  u^  and  u^  of  Eq,  (6,27),  By  virtue 
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of  the  absence  of  a  first  derivative  term  in  the  latter,  F  is  independ¬ 
ent  of  X  and  is  found,  following  McKelvey,^^  to  be  proportional  to 
e"^.  Since  n(x) ,  Eq.  (6.7),  is  pure  imaginary  and  is  a  polynomial 
in  e"^,  the  (n  +  2)-by-(n  +  2)  block  of  Y“^  is  a  bounded  function  of  e. 
Since  r(x,  e),  Eq.  (D.l),  has  a  power  series  expansion  at  l/e  =  0, 
it  is  bounded  in  e  and  the  behavior  of  each  element  of  the  kernel  of 
the  differential  equation  is  known. 

One  now  can  apply  the  theorem  (D.2)  on  any  finite  x  interval, 
say  (O,  2x:*),  which  includes  the  turning  point.  Since  the  behavior  of 
the  y^  and  y„  columns  of  Y  depends  critically  on  the  value  of  x,  we 
divide  this  interval  into  three  subintervals  (O, 

(x_^,  2x^)  by  the  condition 

lt[z(xj,  ie]  1  =  M  (D.4) 

with  M  some  large  positive  number.  Since  f  is  linear  in  e  and  of  order 

(x  -  x^)  near  the  turning  point,  the  central  subinterval,  containing 

1 

the  point  x*,  has  length  of  order  e 

By  virtue  of  this  decomposition,  we  can  apply  the  z  >  0  asymptotic 
formulas  for  y^  and  Eq.  (6.iK)a)  and  (6.4la),  on  the  interval 

(x^,  2x*) ,  and  the  z  <  0  formulas  Eq,  (6,40b)  and  (6.4lb)  on  (O,  x_) . 

On  the  central  subinterval,  however,  the  asymptotic  formulas  no  longer 
apply  so  that  we  return  to  the  full  expression  for  the  u^(z,  co) functions 
in  Eq,  (6,29).  Since  f  is  bounded,  it  follows  from  the  power  series 
Eq.  (6.31)  that  the  v^(z,  £)  functions  in  Eq.  (6.29e)  are  also  bounded. 
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In  addition,  it  can  be  shown  from  this  power  series  that,  although  the 
derivatives  of  the  confluent  hypergeometric  functions  are  infinite  at 
f  =  0  (i.e.,  at  the  turning  point),  it  is  nonetheless  true  that  v^ 
and  vi^  are  bounded  throughout  the  subinterval  x^)  .  Thus  u^  and  its 

derivative,  which  determine  y^  and  through  Eq.  (6.59) j  are  also 
bounded . 

To  prove  admissibility  of  either  y^  or  we  begin  on  the 
subinterval  on  which  it  is  recessive,  setting  the  lower  integration  limit, 
for  all  matrix  elements,  in  Eq.  (D.l)  to  the  outer  limit  of  the  interval, 
i.e.,  either  0  or  2x*.  For  the  sake  of  illustration,  we  will  consider 
the  admissibility  of  y^,  so  that  we  begin  on  (x_^,  2x^) .  We  set  the 
function  h(x,  e)  of  the  theorem  to  ^  so  that,  by  Eq.  (6.40a),  the 
order  relation  for  the  inhomogeneous  part  in  Eq.  (D.2)  holds.  The 
required  order  relation  for  the  kernel  can  then  be  established  through 
the  asymptotic  expressions  for  and  y^  and  the  boundedness  of  y^^  y|^3 
thereby  proving  admissibility  on  one  subinterval. 

On  the  central  interval,  the  function  h(x,  e)  in  the  theorem 
becomes  simply  unity,  and  the  integral  is  split  into  a  contribution 
over  the  outer  subinterval  (x_^,  2x^)  and  one  from  x^  to  x.  The  first 
contribution  involves  only  known  functions,  since  on  2x*)  was 

just  found  to  be  y^  in  the  last  paragraph,  and  is  grouped  with  the 
inhomogeneous  part  of  Eq.  (D.l)  to  form  the  inhomogeneous  part  g(x,  e) 
of  the  theorem.  The  order  relations  are  once  more  established  using 
the  boundedness  of  the  y^  within  the  (x  ,  x_^)  and  the  fact  that  the 
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interval  of  integration  is  of  order  e"^.  Thus  is  proved  equal  to 
to  order  e“”^,  on  (x^,  2x*) , 

The  final  interval  (0,  x  )  is  considered  by  combining  the  integral 

from  X  to  2x-?«-  in  Eq.  (D.l)  with  y  to  fom  g(x,  e)  for  Eq.  (0.2).  One 

^  -k-- 

shows  that  this  part  satisfies  the  first  of  Eq.  (D,2b)  for  h  =  ^ 

and  then  turns  to  the  inhomogeneous  part,  using  again  the  known 

a,symptotic  e3<pressions.  The  proof  of  admissibility  is  then  complete. 
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